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ABSTRACT 

Complete sequences of new analytic solutions of Einstein's equations which describe 
thin super massive disks are constructed. These solutions are derived geometrically. The 
identification of points across two symmetrical cuts through a vacuum solution of Einstein's 
equations defines the gradient discontinuity from which the properties of the disk can be 
deduced. The subset of possible cuts which lead to physical solutions is presented. At 
large distances, all these disks become Newtonian, but in their central regions they exhibit 
relativistic features such as velocities close that of light, and large redshifts. Sections with 
zero extrinsic curvature yield cold disks. Curved sections may induce disks which are stable 
against radial instability The general counter rotating flat disk with planar pressure tensor 
is found. Owing to gravomagnetic forces, there is no systematic method of constructing 
vacuum stationary fields for which the non-diagonal component of the metric is a free 
parameter. However, all static vacuum solutions may be extended to fully stationary 
fields via simple algebraic transformations. Such disks can generate a great variety of 
different metrics including Kerr's metric with any ratio of a to m. A simple inversion 
formula is given which yields all distribution functions compatible with the characteristics 
of the flow, providing formally a complete description of the stellar dynamics of flattened 
relativistic disks. It is illustrated for the Kerr disk. 
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1. Introduction 

The general-relativistic theory of rapidly rotating objects is of great intrinsic interest and has poten- 
tially important applications in astrophysics. Following Einstein's early work on relativistic collisionless 
spherical shells, Fackerell (1968) g§, Ipser (1969) §J| and Ipser & Thorne (1968) f| have developed the 
general theory of relativistic collisionless spherical equilibria. The stability of such bodies has been further 
expanded by Katz (1975) jy]]. Here the complete dynamics of flat disk equilibria is developed and il- 
lustrated, extending these results to differentially rotating configurations with partial pressure support. A 
general inversion method for the corresponding distribution functions is presented, yielding a coherent model 
for the stellar dynamics of these disks. Differentially rotating flat disks in which centrifugal force almost 
balances gravity can also give rise to relatively long lived configurations with large binding energies. Such 
objects may therefore correspond to possible models for the latest stage of the collapse of a proto-quasar. 
Analytic calculations of their structure have been carried into the post Newtonian regime by Chandrasekhar. 
Strongly relativistic bodies have been studied numerically following the pioneering work on uniformly ro- 
tating cold disks by Bardeen & Wagoner [fij. An analytic vacuum solution to the Einstein equations, the 
Kerr metric fll2|] , which is asymptotically flat and has the general properties expected of an exterior metric 
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of a rotating object, has been known for thirty years, but attempts to fit an interior solution to its exterior 
metric have been unsatisfactory. A method for generating families of self-gravitating rapidly rotating disks 
purely by geometrical methods is presented. Known vacuum solutions of Einstein's equations are used to 
produce the corresponding relativistic disks. One family includes interior solutions for the Kerr metric. 

1.1. Background 

In 1919 Weyl (2(| and Levi-Civita gave a method for finding all solutions of the axially symmetric Einstein 
field equations after imposing the simplifying conditions that they describe a static vacuum field. Bicak, 
Lynden-Bell and Kalnajs (1993) - here-after refered as BLK - derived complete solutions corresponding 
to counter rotating pressure-less axisymmetric disks by matching exterior solutions of the Weyl- Levi-Civita 
type on each side of the disk. For static metrics, the vacuum solution was obtained via line superposition 
of fictitious sources placed symmetrically on each side of the disk by analogy with the classical method of 
images. The corresponding solutions describe two counter-rotating stellar streams. This method will here be 
generalised, allowing true rotation and radial pressure support for solutions describing self-gravitating disks. 
These solutions provide physical sources (interior solutions) for stationary axisymmetric gravitational fields. 

1.1.1. Pressure Support via Curvature 

In Newtonian theory, once a solution for the vacuum field of Laplace's equation is known, it is straight- 
forward to build solutions to Poisson's equation for disks by using the method of mirror images. In general 
relativity, the overall picture is similar. In fact, it was pointed out by Morgan & Morgan (1969) Jl^] that for 
pressureless counter rotating disks in Weyl's co-ordinates, the Einstein equation reduces essentially to solving 
Laplace's equation. When pressure is present, no global set of Weyl co-ordinates exists. The condition for 
the existence of such a set of co-ordinates is that Trr + T zz = 0. On the disk itself, this is satisfied only 
if the radial eigenvalue of the pressure tensor, pr, vanishes (T zz vanishes provided the disk is thin). When 
Pr is non-zero, vacuum co-ordinates of Weyl's type still apply both above and below the disk, but as two 
separate co-ordinate patches. Consider global axisymmetric co-ordinates (R',cj)',z') with z' = on the disk 
itself. For z' > 0, Weyl's co-ordinates R, 4>, z may be used. In terms of these variables, the upper boundary 
of the disk z' = + is mapped onto a given surface z — f(R). On the lower patch, by symmetry, the disk will 
be located on the surface z = —f(R). Thus in Weyl co-ordinates, the points z = f(R) on the upper surface 
of the disk must be identified with the points z = —f(R) on the lower surface. The intrinsic curvature of the 
two surfaces that are identified match by symmetry. The jump in the extrinsic curvature gives the surface 
distribution of stress-energy on the disk. A given metric of Weyl's form is a solution of the empty space 
Einstein equations, and does not give a complete specification of its sources (for example, a Schwarzschild 
metric of any given mass can be generated by a static spherical shell of any radius). With this method, 
all physical properties of the source are entirely characterised once it is specified that the source lie on the 
surface z = ±f(R) and the corresponding extrinsic curvature (i.e. its relative layout within the embedding 
vacuum spacetime) is known. The freedom of choice of v in the Weyl metric corresponds to the freedom 
to choose the density of the disk as a function of radius. The supplementary degree of freedom involved in 
choosing the surface of section z = f(R) corresponds classically to the choice of the radial pressure profile. 

1.1.2. Rotating Disks 

For rotating disks, the dragging of inertial frames induces strongly non-linear fields outside the disk which 
prohibit the construction of a vacuum solution by superposition of the line densities of fictitious sources, 
nevertheless, in practise quite a few vacuum stationary solutions have been given in the literature, the most 
famous being the Kerr metric. Hoenselaers, Kinnersley & Xanthopoulos (1979) jjl|[l4]] ( HKX hereafter) 
have also given a discrete method of generating rotating solutions from known static Weyl solutions. The 
disks are derived by taking a cut through such a field above all singularities or sources and identifying it 
with a symmetrical cut through a reflection of the source field. This method applies directly to the non- 
linear fields such as those generated by a rotating metric of the Weyl-Papapetrou type describing stationary 
axisymmetric vacuum solutions. The analogy with electromagnetism is then to consider the field associated 
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with a known azimuthal vector potential A^, the analog of gt^/gtt- The jump in the tangential component 
of the corresponding .B-field gives the electric current in the plane, just as the jump in the t, <j) component 
of the extrinsic curvature gives the matter current. 

The procedure is then the following: for a given surface embedded in curved space-time and a given 
field outside the disk, the jump in the extrinsic curvature on each side of the mirror images of that surface 
is determined in order to specify the matter distribution within the disk. The formal relationship between 
extrinsic curvatures and surface energy tensor per unit area was introduced by Israel and is described in 
details by Misner Thorne Wheeler (1973) p 552. Plane surfaces and non-rotating vacuum fields (ie.the 
direct counterpart of the classical case) were investigated by BLK and led to pressureless counter-rotating 
disks. Any curved surface will therefore produce disks with some radial pressure support. Any vacuum 
metric with g t( f, ^ outside the disk will induce rotating disks. 



Figure 1.1: Symbolic representation of a section z = ±/(i?) through the gttj, field 
representing the upper and lower surface of the disk embedded in two Weyl-Papapetrou 
fields. The isocontours of lu(R, z) correspond to a measure of the amount of rotation in 
the disk. 



In section 2, the construction scheme for rotating disks with pressure support is presented. In section 3, 
the properties of these disks are analysed, and a reasonable cut, z = f(R), is suggested. Section 4 describes 
stellar equilibria which have the corresponding anisotropic stress energy tensors, while section 5 derives the 
properties of the limiting counter-rotating disks. This method is then first applied to warm counter rotating 
disks in section 6; the dominant energy condition given for Curzon disks by Chamorow et. al. (1987) [EJ, and 
Lemos ' solution for the self - similar Mestel disk with pressure support are recovered. In section 7, the Kerr 
disk is studied and the method for producing general HKX disks is sketched. 



2. Derivation 

The jump in extrinsic curvature of a given profile is calculated and related to the stress energy tensor 
of the matter distribution in the corresponding self gravitating relativistic disk. 



2.1. The Extrinsic Curvature 

The line element corresponding to an axisymmetric stationary vacuum gravitational field is given by 
the Weyl-Papapetrou (WP) metric 



ds 2 



e 2 v (dt - u d<\>) 2 + e 2 c ~ 2 v {dR 2 + dz 2 ) + R 2 



-1v 



(2.1) 



where (R, <ft, z) are standard cylindrical co-ordinates, and v, £ and u> are functions of (R, z) onlygeo- 
metric units G = c = 1 are used throughout this paper. Note that this form is explicitly symmetric under 
simultaneous change of <f> and t. The vacuum field induces a natural metric on a 3-space-like z = f(R) 
surface 

da 2 = q a g — —rdx a dx b = q a g h"h? dx a dx b , 

p dx a dx b p a b 



(2.2) 



7 ab dx a dx b 



where {x a } are the co-ordinates on the embedded hypersurface z = f(R), {x a } = (t,R,(f>), and {x a } = 
(t, R, z, 4>) are Weyl's co-ordinates for the embedding spacetime. Here, jab stands for the embedded metric, 
and h" is given by 

"10 

10/', (2.3) 
10 



dx a 
dx a 
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where hereafter a prime, ( )', represents the derivative with respect to R. The line element is then 

da 2 =-e 2 » (dt~wd<pf + R 2 e- 2v dcjy 2 + e 2 ^ 2 » (l + J' 2 ) dR 2 . (2.4) 
Let the upward pointing normal AT„ to the surface z = f(R) be 

d 



N » = ^r[*-/(fl)] = «W'.o,i). 



The normalised normal is therefore 



N u 



N, 



The extrinsic curvature tensor K is given in its covariant form by 



(d a 



1 a/3 a a n b I > 



(2.5) 
(2.6) 

(2.7) 



K tt = 

K RR 

= 

Ka,a = 



dv 



where i9 a /i^ = dh^/dx a . Eq. (2.6)| and (2.3) together with the Christoffel symbols for the Weyl metric lead 
to the computation of K , the extrinsic curvature of that surface embedded in the Weyl Papapetrou metric: 

e-< {9u\ (2ga) 

(2.86) 
(2.8c) 
(2.8d) 



Vi + / 



f" 



0("-O 

dN 



) (i + /' 2 ) , 



d[u+i log (w) 



ON 
R dN) 



uie 



where the notation d/dN = (djdz — f'd/dR) has been used. 



log (w)] 



dN 



2.2. The Stress Energy Tensor 



The stress energy tensor per unit surface t£ is the integral of of the stress energy tensor carried along 
the normal to the surface z = f(R). In a locally Minkowskian frame co- moving with the mean flow of the 
disk, the corresponding orthonormal tetrad is 



er(l,0,-w). 



=w = e f-Vi + /' 2 (o,i,o), 

Re'" (0,0,1), 



-,(2) 



so that 



ds 2 = 77(a)(6) (efW) (efW) 
with 77(a)(6) = Diag(— 1, 1, 1). In that frame, 



»(6) 



£ 

p R 
vs 



(2.9a) 
(2.96) 
(2.9c) 

(2.10) 



(2.11) 



After a Lorentz transformation to a more general frame in which the flow is rotating with relative velocity 
V in the 6 direction, this stress becomes 



»(f>) _ 



1 



1 - V 2 



e+p^V 2 

(1 - V 2 ) p R 

(p + e) V 



, + e) V 


+ eV 2 



(2.12) 
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2.3. The Discontinuity Equations 

Israel (1964) has shown that Einstein's equations integrated through a given surface of discontinuity 
can be re-written in terms of the jump in extrinsic curvature, namely 



T a — — \K a — K a ^ a i + = r a 



(2.13) 



where [ ]_ stands for ( ) taken on z = f(R) minus ( ) taken on z — —f(R). This tensor C% is known as 
the Lanczos tensor. 

In the tetrad frame (2.9), the Lanczos tensor given by Eqs. (2.13) and (2.8) reads 



c (a)(b) = 



47r v /IT7 7i 



r T dN 



au) e 

dN 2R 





_£ 

R 




ou e 

dN 2R 



i+/' 2 ~ aw 



(2.14) 



Identifying the stress energy tensor with the tetrad Lanczos tensor £( a '( b ) according to Eq. (2.13), 

and solving for p R ,p#,e and V yields 



V 



Re 



-2i/ 



duj/dN 



r 



4^ v /i+7 



4 7 r v /TT7 7 2 



4^1 + /' 2 V ^ 







+ ( r 


f 


\i + r 2 


+ 2R + - 


( r 


+ ^ + 
2R 


\i + r 2 


-r\ 




r r 





dN 



dN 



(2.15a) 
(2.156) 
(2.15c) 
(2.15d) 



where 



dv \ 2 e 4 " / dtu x J 



R dN R 2 V dN 



(2.16) 



All quantities are to be evaluated along z = f(R). Eq. (2.15) gives the form of the most general solution to 
the relativistic rotating thin disk problem provided the expressions for e,p^,p R are physically acceptable. 

3. Physical properties of the warm disks 

Physical properties of interest for these disks are derived and related to the choice of profile z = f(R) 
compatible with their dynamical stability. 

Defining the circumferential radius R c , proper radial length R, and a synchronised proper time t* by 



R r = Re'" 



R = J + f' 2 e c - u dR, 



(3.1a) 
(3.16) 

(3.1c) 



where the integral over dR is performed at z = f{R) and that over dt is performed along a given trajectory 
(T), the line element on the disk (2.4) reads 



da 1 



It 2 + dR 2 + R 2 d6 2 . 



(3.2) 
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For circular flows (j^t* = e"[l — u>d<j)/dt] dt), V may be re- written as 

V = R C ^. (3.3) 

Equation (3.3) is inverted to yield the angular velocity of the flow as measured at infinity: 

# = Ve"/R c 

dt l + uVe»/R c - [ ' 

Similarly, the covariant specific angular momentum h of a given particle reads in terms of these variables 

h = ^= lttf> u t + 1 ^u* = ^=L=(R c V + uen , (3.5) 

where /i is the rest mass of the particle. The covariant specific energy e of that particle reads 

e= -- = - ltd> u^-j tt u t = / e " . (3.6) 

n ' ' vi - v 2 

which gives the central redshift, 1 + z c — exp(— v). The velocity of zero angular momentum observers (so 
called ZAMOs) follows from [Eg. (3.5)[ 

V. = -*£. (3-7) 

With these observers, the cuts z = f(R) may be extended inside ergoregions where the dragging of inertial 
frames induces apparent superluminous motions as measured by locally static observers. The circumferential 
radius R z as measured by ZAMOs is Lorentz contracted with respect to R c , becoming 

Rz = = Rc[l- V?] 1/2 , (3.8) 
while the velocity flow measured by ZAMOs is 

The total angular momentum of the disk H follows from the asymptotic behaviour of the vacuum field uj 
at large radii, to — > —2H/r. Define the binding energy of the disk AE as the difference between the baryonic 
rest mass Mq and the total mass-energy of the disk M as measured from infinity given by v — > —M/r: 

AE = M - M. (3.10) 

M is also given by the Tolman formula, but corresponds, by construction, to the mass of the line source for 
the vacuum field. The baryonic rest mass can be expressed as 

M = J^^gU^dS^ = J ^J^ya ZlirRcdR, (3.11) 

where the baryonic energy density E can be related to the energy density e through some yet unspecified 
equation of state or the knowledge of a distribution function. 
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3.1. Constraints on the internal solutions 

The choice of the profile z = f(R) is open, provided that it leads to meaningful physical quantities. 
Indeed it must lead to solutions which satisfy e > 0, p R >Q,p<p> 0, andl^ < 1 while the dominant energy 
condition implies also that e > p,p and e > p R . These translate into the following constraints on / 



e>0 => ^Af 2 - O 2 + (Af - 2Z) > ; (3.12a) 

0<V < 1 ^7V>OandO>0; (3.126) 
O<P0<£ => AZ Af > (O 2 + AZ 2 ) andW>2Z; (3.12c) 

< p R < e => /' < and - 2f'/R < ^Af 2 - O 2 + (AT - 2Z) ; (3.12d) 

given 

The existence of a solution satisfying this set of constraints can be demonstrated as follows: in the limit 
of zero pressure and counter-rotation (i.e. O — > 0, Af — » 2dv/dz, and Z — > d(/dz), any cut z = Const = 
6 ^> m satisfies Eqs. (3.12). By continuity, there exist solutions with proper rotation and partial pressure 
support. In practise, all solutions given in sections 6 and 7 satisfy the constraints (3.12). Note that in the 
limit of zero radial pressure, Eq. (3.126) implies d/dz [co/R + exp2^] < 0. 

3.2. Ansatz for the profile 

In the following discussion, the section z = f(R) is chosen so that the corresponding radial pressure gradient 
balances a given fraction of the gravitational field which would have occurred had there been no radial 
pressure within the disk. This choice 'bootstraps' calculations for all relevant physical quantities in terms of 
a single degree of freedom (i.e. this fraction rj), rather than a complete functional. This gives 

On the r.h.s. of Eq. (3.14), p a is put to zero and V,p0,e are re-expressed in terms of (, v, u via Eqs. (2.15) 
with /' = /" = and z = Const. On the l.h.s. of Eq. (3.14), p R is chosen according to Eq. (2.15d). In 
practise, this Ansatz for / is a convenient way to investigate a parameter space which is likely to be stable 
with respect to ring formation as discussed in the next subsection. In principle, a cut, /, could be chosen so 
as to provide a closed bounded symmetrical surface containing all fictitious sources. Indeed, by symmetry 
no flux then crosses the z = plane beyond the cut, and the energy distribution is therefore bound to the 
edge of that surface: this corresponds to a finite pressure supported disk. It turns out that in general this 
choice is not compatible with all the constraints enumerated in the previous section. More specifically, the 
positivity of fails for all finite disk models constructed. The Ansatz described above gives for instance an 
upper bound on the height of that surface when assuming that all the support is provided by radial pressure. 
This height is in turn not compatible with positive azimuthal pressure at all radii. 
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3.3. Stability 

To what extent are the equilibria studied in the previous sections likely to be stable under the action of 
disturbances? The basic instabilities can be categorised as follows: 

dynamical instabilities, which are intrinsic to the dynamical parameters of the disks, grow on an orbital 
time scale, and typically have drastic effects on the structure of the system. 

secular instabilities, which arise owing to dissipative mechanisms such as viscosity or gravitational 
radiation, grow on a time scale which depends on the strength of the dissipative mechanism involved, and 
slowly drive the system along a sequence of dynamical equilibria. 

Amongst dynamical instabilities, kinematical instabilities correspond to the instability of circular orbits 
to small perturbations, and collective instabilities arises from the formation of growing waves triggered by 
the self-gravity of the disk. Rings, for instance, will be generated spontaneously in the disk if the local radial 
pressure is insufficient to counteract the self-gravity of small density enhancements. Even dynamically stable 
non-axisymmctric modes may drive the system away from its equilibrium by radiation of gravitational waves 
which will slowly remove angular momentum from the disk. For gaseous disks, viscosity and photon pressure 
will affect the equilibria. Radiative emission may disrupt or broaden the disk if the radiation pressure 
exceeds the Eddington limit. The energy loss by viscosity may induce a radial flow in the disk. However, 
the disks discussed in this paper have anisotropic pressures inappropriate for gaseous disks and the accurate 
description of the latter two processes requires some prescription for the dissipative processes in the gas. 
The scope of this section is therefore restricted to a simple analysis of the dynamical instabilities. 

Turning briefly to the corresponding Newtonian problem, Toomre (1963) |p4| gave the local criterion 
for radial collective instability of stellar disks, 

a B > , (3.15) 

K 

where So is the local surface density, a 2 R the radial velocity dispersion, and k the epicyclic frequency of the 
unperturbed stars. This criterion is derived from the first critical growing mode of the dispersion relation 
for radial waves. The corresponding critical wavelength is Jeans length ~ 2<7^/GXo- 

For the relativistic disks described in this paper, spacetime is locally flat, which suggests a direct 
translation of Eq. (3.15) term by term. The constraints that stability against ring formation places on these 
models will then be addressed at least qualitatively via the Newtonian approach. The proper relativistic 
analysis is left to further investigation. The relativistic surface density generalising Eo is taken to be the 
co-moving energy density s given by Eq. (2.156). The radial velocity dispersion a 2 is approximated by p R /s- 
The epicyclic frequency is calculated in the appendix. Putting Eqs. (DA) and (D.5) into Eq. (3.15) give 
another constraint on / for local radial stability of these disks. Note that the kinematical stability of circular 
orbits follows from Eq. (D.3) by requiring k 2 to be positive. 

4. Stellar dynamical equilibria for rotating super-massive disks 

The method described in section 2 will generally induce rotating disks with anisotropic pressures (p^ ^ Pr)- 
These objects may therefore be described in the context of stellar dynamics. It should be then checked that 
there exist a stellar equilibria compatible with a given vacuum field and a given cut z = f(R). 
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4.1. Distribution functions for rotating super-massive disks 

Here a general procedure to derive all distribution functions corresponding to a specified stress energy tensor 
is presented for disks with non-zero mean rotation (a more direct derivation for counter-rotating disk is given 
in Appendix B). The detailed description of the dynamics of the disk follows. 
In Vierbein components, the stress energy tensor reads 



T (aj3) _ 



U (e, h) P {a) P {l3) 



(4.1) 



where /+ (e, h) is the distribution of stars at position R, <j> with momentum p( R > ,p(w. For a stationary disk, 
it is a function of the invariant of the motion e, h. Now for the line element (2.4): 



do 2 = - e 2v {dt - uj d(j)f + e 2 ^ 2 " (l + f 2 ) dR 2 + R 



> 2 e- 2v d6 2 , 



the Vierbein momenta read 



=e v R- 1 (h-ue) 
PW = e -"e, 



p(R) 



e 2g-2i/ „2u d-2 



e 2u R~ 2 (fc-we) -1 



1/2 



Calling x = h/e and i? = 1/e 2 , Eqs. (4.3) becomes 



p(K) = 0-1/2 



,-2j/ „2v d-2 



1/2 



In terms of these new variables, the integral element dP^dP^ becomes 



dP W dP (R) 



gp(<P) dP^ 



dxdd = 



e- v R- Y $- 2 dxdd 



2\je- 2v - e 2v R- 2 (x -(*>) -# 
Given Eq. (4.5)| and Eq. (4.36), Eq. (4.1) may be rewritten as 



rM)= f fpMpiff) h{e,h)R-^-*/ 2 dxdi) 



In particular, 



RT (tt) e 2v 



p2 Ji(t4>) _ 



p3 j^((f>4>) e ~2v 



2\Je~ 2v - e 2v R- 2 (x - w) - d 

A(e,/Q $- 5 / 2 d X d'& 
2 \/ e- 2u - e 2v R- 2 ( X - oj) 2 - ■& 
U(e,h) iT 5 / 2 dxdtf 



(x-w) 



2yje- 2v - e 2v R- 2 { X - w) - 
, U{e,h) $-">l 2 dxd§ 



2\je- 2v - e 2 " R- 2 { X - w) - 
Note that given Eqs. (4.7), T^ RR > follows from the equation of radial support. Note also that 

U(e,h) e- u tf- 3 / 2 d X d§ 



R, 



2a e 



-2v _ e 2v R -2 ( x _ w ) z _^ 



(4.2) 



(4.3a) 
(4.36) 

(4.3c) 



(4.4a) 
(4.46) 

(4.4c) 



(4.5) 



(4.6) 



(4.7a) 
(4.76) 
(4.7c) 



(4.8) 
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or equivalently, 

T (tt) _ T (00) _ T (RR) _ J J f* ( e M dp W dp (R) _ (- 4 9 ) 

Equation (4.9) yields the total gravitating mass. It is known as Tolman's formulae and can be derived 

directly from the relation P M P M = —1 and Eq. (4.1). 

Let 

f<x,b>-\[ . -fig&r- <«»' 

2J yje-*" - e 2 " i?- 2 ( X - iof - { VJ 7 - ^ 

where /* (x, i?) = /* (e, h) iT 5/2 , and y (R, X ) = - e 2 » R~ 2 ( X ~ ujf . (4.11) 
The set of Eqs. (4.7) becomes 

RTW e 2v = j F ( x , R) d X = Re~ 2v J F (v+,R) dvj, , (4.12a) 

R?T {t ® = f F( X ,R) (x-u)dx =R 2 e- 4v J v^F{v^R) dv,p , (4.126) 
R 3 e~ 2v = J Fix, R) (X - d X = R 3 e- 6 » Jv 2 F ( V4> , R) dv^ . (4.12c) 

The r.h.s of Eqs. (4.12) was re-expressed convinienlty in terms of the new variable = e v (x— u>)/R c , and 
the new function F{v^,, R) = F(x — R c v < f,e~ I/ +u, R) in order to stress the analogy with the classical identities. 
Note that in the limit of zero pressures, reduces to V, the geodesic velocity of stars on circular orbits. 
In Eq. (4.12), the dragging of incrtial frames requires to fix simultaneously three moments of the velocity 
distribution to account for the given energy density, pressures and rotation law. This third constraint does 
not hold in the corresponding Newtonian problem, nor for static relativistic disks as described in appendix B. 
Any positive function F(x, R) satisfying the moment constraints Eqs. (4.12) corresponds to a possible choice. 
For instance the parametrisation is carried in Appendix C for F distributions corresponding to powers of 
Lorentzians in azimuthal momenta. Let F(x, y) — F(x,R), where y(R, x) is given by Eq. (4.11). Written 



in terms of F, the integral equation Eq. (4.10) is solved for /* by an Abel transform 



l 

The latter integration may be carried in R space and yields 

dR 



r _ . . dF\ dR 2 f fdF\ dR lA 

/4e,/ l ) = _/ _ - — — — / — (4.14) 



'R, 



where p( fl ) is given by Eq. (4.3c) as a function of R, e, and h. Here F is chosen to satisfy Eq. (4.12) which 
specifies the stress energy components T^ a ^ . The integration limits correspond to two branches: the first 
branch is to be carried between lZ e (h), the inner radius of a star on a "parabolic" orbit with momentum 
h, and lZ p (h,e), the perigee of a star with invariants {h, e); the second branch contributes negatively to 
Eq. (4.14) and corresponds to radii larger than lZ a (h,e), the apogee of a star with invariants (h,e). Note 
that the derivative in Eq. (4.14) is taken at constant reduced momentum X = h/e. This inversion procedure 
is illustrated on Fig. 7.5| . All properties of the flow follow in turn from /*(e, h). For instance, the rest mass 



surface density, E, may be evaluated from the detailed "microscopic" behaviour of the stars, leading to an 
estimate of the binding energy of the stellar cluster. 
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Figure 4.1: sketch for the effective potential Y(R) as a function of radius. The horizontal 
line corresponds to the energy level 1/e 2 , the absisse line corresponds to the escape energy 
1/e 2 = 1. The ordinate axis corresponds to TZ e (h,e), the inner radius of a star on a 
"parabolic" orbit with momentum h, the left vertical line corresponds to lZ p (h,e), the 
perigee of a star with invariants (h, e); the right vertical line corresponds to lZ a (h,e), 
the apogee of a star with invariants (h,e). The area between the curve Y = Y(R) and 
the lines Y = 1/e 2 and R — TZ e are shaded, defining three regions from left to right. 
The middle region does not contribute to Eq. (4.14). The equation Y — 1/e 2 has two 
roots corresponding to lZ p (h, e) and lZ a (h, e), while Y = 1 has two roots corresponding to 
infinity and lZ e . The sign of the slope of Y(R) gives the sign of the contribution for each 
branch in Eq. (4.14). 



4.2. A simple equation of state 

An equation of state for these rotating disks with planar anisotropic pressure tensors is alternatively found 
directly while assuming (somewhat arbitrarily) that the fluid corresponds to the superposition of two isotropic 
flows going in opposite directions. In other words, the anisotropy of the pressures measured in the frame 
co-moving with the mean flow V is itself induced by the counter rotation of two isotropic streams. (These 
counter rotating streams are described in more details in the next section.) In each stream, it is assumed 
that the pressure is isotropic and that the energy is exchanged adiabatically between each volume element. 
The detailed derivation of the corresponding relationship between the surface density and the pressure of 
each stream is given in Appendix A. The above set of assumption yields the following prescription for £ : 



- V4> ~ Pr + V(e ~ P<t> - Pr) ( £ - P<t> + 3 Ph) , 4 ^ 



2 



where Vq — \J (p^ — p R ) / (e + p R ) is the counter rotating velocity measured in the frame co-moving with V 
which induces p^ ^ p R . In the classical limit, £ — * (e —p^)/^/\ — Vq. The binding energy of these rotating 



disks is computed from Eq. (4.15) together with Eqs. (3.11) and (2.15). 



5. Application: warm counter-rotating disks 

For simplicity, warm counter-rotating solutions are presented first, avoiding the non-linearities induced 
by the dragging of inertial frames. These solutions generalise those of BLK, while implementing partial 
pressure support within the disk. Formally this is achieved by putting w(i?, z) identically to zero in Eq. (4.2); 
the metric for the axisymmetric static vacuum solutions given by Weyl is then recovered: 

ds 2 = -e 2v dt 2 + e 2 ^ 2 - (dR 2 + dz 2 )+R 2 e- 2 »d0 2 , (5.1) 

where the functions ((R, z) and v(R, z) are generally of the form (cf. Eq. (BLK 2.21)) 

W{b)db , (5.2a) 
^/R 2 + (\z\+bf 

W (60 W (b 2 ) Z (b u b 2 ) dbi db 2l (5.26) 

with Z given by 



Here 6i and 62 are the distances of two points below the disk's center and r\ and r 2 are the distances 
measured from these points to a point above the disk. Below the disk, Z is given by reflection with respect 
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to the plane z = 0. When the line density of fictitious sources is of the form W(b) oc b~ m (Kuzmin-Toomre), 
or W(b) oc J( m > (b) (Kalnajs-Mestel), the corresponding functions (,v have explicitly been given in Bicak, 
Lynden-Bell & Pichon (1994) B. For the metric of Eq. (5.1), the Lanczos tensor (2.14) reads 



£(*)(*) 

£(fl)(R) 

£(<t>)M 



o»-C, 



L 

R 

-/' 
R 

-/" , ac 
T+7 75 ' ON 



aw 



(5.4) 



Consider again counter-rotating disks made of two equal streams of stars circulating in opposite direc- 
tions around the disk center. The stress energy tensor is then the sum of the stress energy tensor of each 
stream: 





e 





" 










Pr 






Pd>_ 


2 

i-v 2 



so 



~PoV 2 





1 





^o 2 ) Po 
po 






-eoV 2 



(5.5) 



the pressure po in each stream being chosen to be isotropic in the plane. Subscript ( )o represents quantities 
measured for one stream. Expressions for e,p^ and p R follow from identifying Eqs. (5.5) and (5.4) according 
to Eq. (2.13) . Solving for e , Vo, andpo m Eq. (5.5), given Eq. (5.4) and (2.13)| gives 



£o 



Po 



f" 



4vrVTT7^ L 1 + / 



Sn^iTf^ L R 

\dN 1 + f 
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L 

R 



ON 



I" , d(2v-0 



1 + /'■ 



dN 



(5.6a) 
(5.66) 
(5.6c) 



The above solution provides the most general counter rotating disk model with pressure support. Indeed, 
any physical static disk will be characterised entirely by its vacu um field v a nd its radial pressure profile 
p R , which in turn defines W(b) and f(R) uniquely according to Eqs. (5.66) and (BLK 2.33). The other 

properties of the disk are then readily derived. 

The angular frequency and angular momentum of these disks follow from Eq. (3.5) , (3.4) on putting to 
to zero and V to Vq. The epicyclic frequency k given by Eq. (DA) may be recast as 



2 _ e v dh 2 

which relates closely to the classical expression k 2 — R~ 3 dh 2 /dR. 

The Ansatz given by Eq. (3.14) for z — f(R) becomes, after the substitutions e — ► Eq, V 



(f'/R)' = -v 



dv 
dR 
dv dv 
J] d~Rdl 



dv 
~d~z 



19C 
2dz 



R dR- 1 



(5.7) 

Vo, P4, -> 0, 
(5.8a) 
(5.86) 



where z is to be evaluated at 6. Equation (5.86) follows from Eq. (5.8a) given the (^) component of the 
Einstein equation outside the disk. 

The equation of state for a relativistic isentropic 2D flow of counter-rotating identical particles is derived 
in appendix A while relating the perfect fluid stress energy tensor of the flow to the most probable distribution 
function which maximises the Boltzmann entropy. It reads 



e = 2p 



(5.9) 
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Solving for £ gives 



so 



2p Q + \Je - 2p ^e + 2p 



(5.10) 



which in t he classica l limit gives S — > £o~Po- The binding energy of these counter rotating disks is computed 
here from Eq. (5.10) together with Eqs. (3.11) and (5.6). Alternatively, it could be computed by constructing 
distribution functions using the inversion described in Appendix B. 



6. Examples of warm counter-rotating disks 



6.1. Example 1 : The Self-similar disk 



The symmetry of the self-similar disk allows one to reduce the partial differential equations corresponding 
to Einstein's field equations to an ordinary differential equation with respect to the only free parameter 
9 = arctan(i?/z) (Lemos, 89 [|l7)) . Lemos' solution may be recovered and corrected by the method 
presented in sections 2 and 4. Weyl's metric in spherical polar co-ordinates (p, x, (f>) defined in terms of 
(R, z, <f>) by p — \f{R 2 + z 2 ) and cot(x) = z/R is 



da 2 = -e 2 »dr 2 + p 2 
A self-similarity argument led Lemos to the metric 
ds 2 = \ll 



e~ 2v sin 2 { X )d4> 2 + e 2 ^ (%\d X 



e 2P ~ N d6 2 +e z - N 



+d6' 



where P, N, and Z are given by 



P (0) = log [sin (ft + n)6], 

Am 

N(9) =- log [cos (ft + n) 6»/2] 



Z{9) 



k + n 
8n 2 
(ft + n) ; 



log [cos (ft + n) 9/2} + 2 log (ft + n) . 



(6.1) 

(6.2) 

(6.3a) 
(6.36) 

(6.3c) 



Here the </>'s in both metrics have already been identified since they should in both case vary uniformly in 
the range [0, 2tt[. The rest of the identification between the two metrics follows from the transformation 



X = (n + ft) 
p = r n+k 

T = t, 



with the result that v and £ take the form: 



2i/ = iV + 2nlog (r), 
2C = Z - 2 log (n + ft) . 

This can be written in terms of the new variables given in Eqs. (6.4) as 

4n 2 



C 



(k + n)' 
2n 

(k + n)' 



Mi). 



log 



■ log (p) 



(6.4a) 
(6.46) 
(6.4c) 



(6.5a) 
(6.56) 



(6.6a) 
(6.66) 
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The procedure described earlier may now be applied. In order to preserve the self-similarity of the solution 
and match Lemos' solution, f(R) is chosen so that 



Then 



/' (R) = const. = cot (77) say. 



d f d , d \ _ 1 d 
dN \ dz OR J p sin iff) df) 



(6.7) 
(6.8) 



Putting Eqs. (6.6) and [Eq. (6.8)| into Eq. (2.15) gives 



e 

Pn 



r 

47T 



(!) 



2n(k- 



(k+r, 



cot (77) 



2 " fc , tan (77/2) 



2,r 



(k+ 

- cot (77) 
■tan (77/2) 



(6.9) 



with 77 to be evaluated at (n + k)n/2. The above equations are in accordance with the solution found by 
Lemos by directly integrating the Einstein field equations. 



6.2. Example 2: The Curzon disk 



Figure 6.1: the "Toomre Q number" Q = <r R re/(3.36£o) is plotted here as a function of 
circumferential radius R c for different compactness parameter a = b/M when 77 = 1/25. 



The Curzon disk is characterised by the following pair of Weyl functions: 

v = -a/r, C = -a 2 R 2 / (2 r 4 ) , (6.10) 

with r = ^/R 2 + f 2 , given the dimensionless parameter a = M/b (recall that G = c = 1), all lengths being 
expressed in units of b. This disk is the simplest example of the Kuzmin- Toomre family, where W(b) cx 5(b). 
It also corresponds to the building block of the expansion given in Eqs. (5.2). Equations (5.4) and (5.5) then 
imply 



Figure 6.2: the binding energy over the rest mass AE/Mo for Curzon disks, plotted 
against the compactness parameter a = b/M for different ratios Po/sq fixed by 77 = 
0.01,0.16, •• -,0.7. 



^ + £ + £ {2/ (r 3 - ai? 2 ) + [a (i? 2 - J 2 ) - 2r 3 ] Rf}' 

=£■ ■ (6-11) 

- T ^ + ^[2i? 2 /-(i? 2 -/ 2 )i?/'] 

The weak energy conditions e > 0, e + p R > 0, and e + > which follow are in agreement with 
those found by Chamorow et al(1987) ||. Their solution was derived by direct integration of Einstein's 
static equations using the harmonic properties of the supplementary unknown needed to avoid the patch of 
Weyl metric above and below the disk. The cut z = f(R) corresponds to the imaginary part of the complex 
analytic function, the existence of which follows from the harmonicity of that new function. 



e 

Pn 



Airb^l + f' 2 
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The Ansatz (3.14) leads here to the cut 



f(R) = b + j 



26 to 3 (Ab 2 + 7R 2 ) 



A(b 2 + R 2 ) i05 (6 2 + i? 2 ) 



5/2 



In Fig. 6.1 



(5.10) : (6.11 



the ratio of the binding energy of these Curzon disks (which is derived from Eqs. 
and 



(6.12) 



(3.11), 



12)) over the corresponding rest mass Mq is plotted with respect to the compactness 



parameter b/Mg for different ratios po/so measured at the origin. The relative binding energy decreases in 
the most compact configurations because these contain too many unstable orbits (n 2 < 0). A maximum 
ratio of about 5 % is reached. 



7. Examples of rotating disks and an internal solution for the Kerr metric 

The problem of constructing disks with proper rotation and partial pressure support is more complicated 
but physically more appealing than that of counter rotating solutions. It is now illustrated on an internal 
solution for the Kerr metric, and a generalisation of the solutions constructed by BLK and BLP to warm 
rotating solutions is sketched. 



7.1. A Kerr internal solution 

The functions (is, u>) for the Kerr metric in Weyl-Papapetrou form expressed in terms of spheroidal co- 



ordinates (R — s^J [x 2 =F 1)(1 — y 2 ), z = sxy) are given by 



2C 



and 



— m 2 + s 2 x 2 + a 2 y 2 
(to + s x) 2 + a 2 y 2 

—to 2 + s 2 x 2 + a 2 y 2 
s 2 (x 2 =F y 2 ) 

2 a to (to + s x) (l — y 2 ) 



—m z + s 2 x 2 + a z y 



t2 ,,2 



(7.1) 
(7.2) 

(7.3) 



where s = (±(to 2 — a 2 )) 1 / 2 . Here ± corresponds to the choice of prolate (+), or oblate (— ) spheroidal 
co-ordinates corresponding to the cases a < to and a > to respectively. In this set of co-ordinates, the 
prescription given in section 2 leads to completely algebraic solutions. The normal derivative to the surface 
z = f(R) reads (when a > to) 



where 



(d/dN) = N (x, y) (d/dx) + N (y, x) (d/dy) , 



N (x, y) = (x 2 - 1) [y/s + x (y 2 - l) f'/R] / (x 2 - y 2 ) . 



(7.4) 



Differentiating Eqs. (7.1) through (7.3) with respect to x and y together with Eq. (7.4) leads via Eqs. (2.15) to 
all physical characteristics of the Kerr disk in terms of the Kerr metric parameters to and a, and the function 
z = f(R) which must be chosen so as to provide relevant pressures and energy distributions according to 
Eqs. (3.12). (Bicak & Ledvinka (1993) ^ have constructed independently a cold Kerr solution when a < to.) 
On the axis R = 0, / = b, and /" = -c, while Eqs. (3.12) imply /' = Oat (x = b/s, y = 1). This in turn 
implies 



Mo = (Pr)o = 



47T 




2m 


"(6 + to) 2 -a 2 




(a 2 + {b + m) 2 ^ 


(a 2 + b 2 — to 2 ) 



(7.5a) 
(7.56) 
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Figure 7.1: the azimuthal velocity (a), the radial pressure (b), the ratio of radial over 
azimuthal pressure and the relative fraction of Kinetic Support (d) namely eV 2 /(eV 2 + 
P<t> +Ph) f° r the a/m = 0.5 Kerr disk as a function of circumferential radius R c , for a class 
of solutions with decreasing potential compactness, b/m — 1.1, 1.2 •■■ 1.7, and relative 
radial pressure support rj/m = 1/5 + 2(b/m — l)/5 . 



Figure 7.2: as in |Fig. 7.l| with a/m — 0.95, rj/m = 1/15 
relative kinetic support. 



(b - l)/5. Note the large 



Figure 7.3: the azimuthal velocity (a), the azimuthal pressure (b), the relative fraction 
of Kinetic Support (c), namely eV 2 /(eV 2 and the angular momentum (d) for the 

a/m = 10 zero radial pressure Kerr disk as a function of circumferential radius R c , for a 
class of solutions with decreasing potential compactness b/m = 210, 260 ■ ■ ■ 510. The first 
curve corresponds to a cut which induces negative azimuthal pressure (tensions). 



where ( )o stands for ( ) taken at R = 0. The constraint that all physical quantities should remain positive 
implies 



b > \J m 2 ~ a 2 , 



c < 2m 



(b + m) — a 2 



(a 2 + {b + m) 2 J (a 2 + b 2 - m 2 ) 



(7.6a) 
(7.66) 



Equation (7.6a) is the obvious requirement that the surface of section should not enter the horizon of the 
fictitious source. Similarly, ergoregions will arise when the cut z = f(R) enter the torus 



z \ -' 
m / 



R a 1 m 
m m 2 a 



1 m 
2~a 



(7.7) 



The analysis may b e extended bey ond the e rgoregion via the ZAMO frame. The characteristics of these 

and (3.9). The central redshift 



frames are given by Eq. (3.7), (3 



1 



z r . = e = 



(m + b) 2 + ■ 
b 2 + a 2 — rr 



(7.S 



can be made very large for such models constructed when b — ► \fr 



Fig. 7.1 




Fig. 7.2 


and 


Fig. 7.3 



anisotropy of the planar pressure tensor (p^, ^ p R ). The anisotropy follows from the properties of the u> 
vacuum field which gives rise to the circular velocity curve. Indeed, in the outer part of the disk, the specific 
angular momentum h tends asymptotically to a constant. The corresponding centrifugal force is therefore 
insufficient to counter-balance the field generated by the v function. As this construction scheme generates 
such equilibria, the system compensates by increasing its azimuthal pressure away from isotropy. When 
a > m, and rj = the height of the critical cut corresponding to the last disk with positive pressures 
everywhere scales like a (a 2 — m 2 ) 1 / 2 in the range 2 < a/m < 20. Note that this limit is above that of the 
highly relativistic motions which only occurs when b tends asymptotically to the Kerr horizon (a 2 — m 2 ) 1 / 2 . 
The binding energy which follows Eq. (4.15) reaches values as high as one tenth of the rest mass for the most 
compact configurations. The inversion method described in section 4.1 was carried for the Kerr disks and 



yields distribution functions which characterise completely the equilibria as illustrated on Fig. 7.5. It was 
assumed that the velocity distr ibution o f stars at radius R was a modified squared Lorentzian, as discussed 
in Appendix C and shown on Fig. 7.4 Note that disks with low rotation support present in their outer 
parts two distinct maxima, one of which is counter rotating, in agreement with the assumptions yielding the 
equation of state Eq. (4.15). 
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Figure 7.4: Parametrization of the number of star with velocity at radius R as a 
squared Lorentzian (as discussed in Appendix C) for the Kerr disk. The corresponding 
radius is appended on each curve. The left panel correspond to a = 0.2 while the right 
panel corresponds to a = 0.8. The pressure law is given in the caption of Fig. 7.5| where 
the corresponding distribution function is illustrated. 



Figure 7.5: the isocontours of the distribution function for the Kerr disk as a function 
of the reduced momentum \ = h/e, and the relative eccentricity of the orbit e c — (e ~ 
e x )/(l — e x ) where e x is the energy of circular orbits with reduced momentum \- The 
four panels correspond to a — 0,0.1,0.2,0.3,0.4 from left to right and top to bottom, 
while the pressure law was chosen to correspond to the cut f(R) — —4(4+ 7i? 2 )/(525(l + 
i? 2 ) 5 / 2 ) + (16 + 15i? 2 )/(10(l + i? 2 )) corresponding to a tepid disk. The fi gures are centered 
on the inner parts of the disk to illustrate the relative shift of the maximum number of 
stars towards larger angular momentum for faster rotating disks. The inversion assumed 
a squared Lorentzian for the distribution in angular velocity as discussed in Appendix C 
and illustrated on Fig. 7.4. 



7.2. Other rotating solutions 



The rotating disk models given in section 2 require prior knowledge of the corresponding vacuum solutions. 
However, while studying the symmetry group of the stationary axially symmetric Einstein Maxwell equations, 
Hoenselaers, Kinnersley & Xanthopoulos (1979) Jl3|,[l4| found a method of generating systematically complete 
families of stationary Weyl-Papapetrou vacuum metrics from known static Weyl solutions. Besides, the 
decomposition of the Weyl potential v$ into line densities provides a direct and compact method of finding 
solutions to the static field. The N+l rank zero HKX transform is defined as follows: let vq and £o be 
the seed Weyl functions given by Eqs. (5.2). HKX define the N x N matrix T parameterised by N twist 
parameters a k and N poles, ak, k — 1, • • • N as 



(r±) 



k,k' 



ia>k ■ 



Tk - r k > 
a k - aw 



± 1 



(7.9) 



with rfc = r(ctfc) = \J R 2 + {z — ak) 2 and (5k = P(a-k) , P being a function satisfying (7.11a) below. Defining 
r = Diag(rfc), they introduce the auxiliary functions 



D 



and L± = 2D ± tr 



w± = D + ± e^°D~, and M± 
which lead to the Weyl-Papapetrou potentials 



+ r±) 



,2" 



3? 



D-D' 



e 2 ? =/ce 2Co 3£ 

Ml (uj+ + uj-) - M- (v+ + u-Y 
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(7.10a) 
(7.106) 

(7.10c) 



where and 3 stand respectively for the real part and the imaginary part of the argument and ( )* 
represents the complex conjugate of ( ). A: is a constant of integration which is fixed by the boundary 
conditions at infinity. The two yet unspecified function (3 and w satisfy: 



(a k ) = — ([z- a k ] V + i?e x V) u , 

^k 

Vva = Re^ x Vfo, 



(7.11a) 
(7.116) 
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Figure 7.6: the azimuthal velocity for the a/ny zero radial pressure Yamazaki disk as 
a function of circumferential radius i? c , for a class of solutions with decreasing potential 
compactness b/ny = 3.6, 3.8, ■ • • 5. Ky is the natural unit of distance in these disks as 
defined by Yamazaki's Eq. (Y-2) @. 



where V = (d/dR,d/dz) and is the unit vector in the 4> direction. The linearity of these equations 
suggests again that solutions should be sought in terms of the line density, W(b), characterising i/q, namely 



/3(a fc ) = 



W(b) 



R 2 + (z-b-a k ) 
R 2 + (z- bf 



db, 



W(b) jb-z) 
R 2 + (z~ bf 



■.db. 



(7.12a) 
(7.126) 



It is therefore a matter of algebraic substitution to apply the above procedure and construct all non-linear 
stationary vacuum fields of the Papapetrou form from static Weyl fields. Following the prescription described 
in section 2, the corresponding disk in real rotation is then constructed. The requirement for the physical 
source to be a disk is in effect a less stringent condition on the regularity of the vacuum metric in the 
neighbourhood of its singularity because only the half space which does not contain these singularities is 
physically meaningful. A suitable Ehlers transformation on D-/D+ ensures asymptotic flatness at large 
distance from the source. 

To illustrate this prescription consider the vacuum field given by Yamazaki (1981) P7| . This field 
arises from 2 HKX rank zero transformation on the Zipoy-Voorhess P5| static metric given by 



v = §tfloj 



(x + iy 



Co =J<riog 



2, (H 



{x 2 - y 2 ) ' 



(7.13a) 
(7.136) 



in prolate spheroidal co-ordinates. It corresponds to a uniform line density between ±<5. The functions (3 
and w follow from Eqs. (7.12). Choosing poles a^, k = 1, 2 at the end of the rod ±5, leads to 



P(a ± ) = ±5 log 
vo = 2Sy. 



x l -\ 

[x t y] 2 



(7.14a) 
(7.146) 



Equations (7.10) together with (7.13) and (7.11) characterize completely the three Weyl-P apapetro u metric 
functions v, C, andw; all physical properties of the corresponding disks follow. For instance, Fig. 7.6 gives the 
zero radial pressure velocity curve of the Yamazaki disk spanning from the Schwarzschild (S — 1/2) metric 
when ai = ct2- 

More generally, the extension of this work to the construction of disks with planar isotropic pressures 
should be possible by requiring that Eq. (2.15c) is identically equal to Eq. (2.15d) given Eq. (7.10). 
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8. Conclusion 

The general counter rotating disk with partial pressure support has been presented. The suggested method 
for implementing radial pressure support also applies to the construction of stationary axisymmetric disk 
with rotation but requires prior knowledge of the corresponding vacuum solution. In this manner, a disk-like 
source for the Kerr field has been constructed orbit by orbit. The corresponding distribution reduces to a 
Kcplerian flow in the outer part of the disk and presents strong rclativistic features in the inner regions such 
as azimuthal velocities close to that of light, large central redshift and ergoregions. The disk itself is likely to 
be stable against ring formation and the ratio of its binding energy to rest mass can be as large as 1:10. The 
broad lines of how to construct all rotating disks arising from HKX-transforming the corresponding counter 
rotating model into a fully self-consistent model with proper rotation and partial pressure support has been 
sketched. It should be a simple matter to implement this method with the additional requirement that the 
pressure remains isotropic with a sensible polytropic index. One could then analyse the fate of a sequence 
of gaseous disks of increasing compactness and relate it to the formation and evolution of quasars at high 
redshift. Alternatively, the inversion method described in section 4.1 yields a consistent description of the 
detailed dynamics for all disks in terms of stellar dynamics. 
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A. Equation of state for a relativistic 2D adiabatic flow 

The equation of state of a planar adiabatic flow is constructed here by relating the perfect fluid stress 
energy tensor of the flow to the most probable distribution function which maximises the Boltzmann entropy. 
Synge pl|] gives an extensive derivation of the corresponding equation of state for a 3D flow. These properties 
are local characteristics of the flow. It is therefore assumed that all tensorial quantities introduced in this 
section are expressed in the local Vierbein frame. 

Consider a infinitesimal volume element d-d defined in the neighbourhood of a given event, and assume 
that the particles entering this volume element are subject to molecular chaos. The distribution function F 
of particles is defined so that F(R, P) ddd 2 P gives the number of particles in the volume d centered on R 
with 2-momentum pointing to P within d 2 P. 

The most probable distribution function F* for these particles is then given by that which maximises 
Boltzmann entropy S subject to the constraints imposed by the conservation of the total energy momentum 
and by the conservation of the number of particles within that volume. This entropy reads in terms of the 
distribution F 

S= -d§ [Flog Fd 2 P. (A.l) 



The stress energy tensor corresponds to the instantaneous flux density of energy momentum through the 
elementary volume d§ (here a surface) 

f d 2 P 

T ab = / F*P a P< 3 ^— , (A.2) 
J \£ I 

where the l/\£\ factor accounts for the integration over energy-momentum space to be restricted to the 
pseudo-sphere P a P a = m 2 . Indeed the detailed energy-momentum conservation requires the integration to 
be carried along the volume element 

f 5 (P a P a - to 2 ) d 3 P = J 6 (P° 2 - £ 2 ) d 2 PdP° = ^ , (A3) 



where £ = \J m 2 + P 2 . It is assumed here that momentum space is locally flat. Similarly the numerical flux 
vector reads 

d 2 P 

F*P Q — . (A4) 

The conservation of energy-momentum of the volume clement di) implies that the flux of energy momentum 
across that volume should be conserved; this flux reads 

T af3 rip d$ = d$ J F*P a d 2 P = Const. , (A5) 

where rip is the unit time axis vector. Keeping the population number constant provides the last constraint 
on the possible variations of S 



cf%n a dd = dz9 / F*d z P = Const. . (A6) 

Varying Eq. (A.l) subject to the constraints (A. 5) (A. 6), and putting SS to zero leads to 

(log F* + 1)8F* = a6F* + X a P a 5F* , (A.7) 

where a and A Q are Lagrangian multipliers corresponding resp. to (A. 5) (A. 6). These multipliers are 
independent of P a but in general will be a function of position. The distribution which extremises S 
therefore reads 

F* = Cexp (X a P a ) . (A.8) 
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The 4 constants C and A a are in principle fixed by Eqs. (A5) and (A 6). The requirement for F* to be 
Lorentz invariant - that is independent of the choice of normal n a - and self consistent, is met instead by 
demanding that F* obeys Eq. (A 2) and (A. 4) , namely 



C J P<*exp(X^)^- = €, 

C J P a P exp (X^) ^f- = T a ?. 



where </>" and T a " satisfy in turn the covariant constraints 



dx c 



and 
0. 



(A9) 
(A10) 

(Alia) 
(A116) 



Equations (A9)-(A11) provide a set of thirteen equations to constraint the thirteen functions C, X a , 0* , and T af3 . 
Equations (A9) and (A 10) may be rearranged as 



<9$ R <9 2 $ 

_ _____ rpotp 

~ dX a ' ~ dX a dX? ' 



(A12) 



where the auxiliary function $ is defined as $ = J exp(A Q P a ) d 2 P/\£\; $ is best evaluated using pseudo polar 
coordinates corresponding to the symmetry imposed by the energy momentum conservation P a P a — m 2 = 0; 
these are 

P 1 = to sinh \ cos <f> , 

P 2 = to sinh % sin , (A. 13) 

P° = im cosh x , 

where the coordinate x is chosen so that the normal n a lies along \ = 0. In terms of these variables, d 2 P/\£\ 
then reads 

d 2 P/\£\ = dP 1 dP 2 /mcoshx = msinhxdx#- (A14) 

Moving to a temporary frame in which A is the time axis ( A = iX — i(— AqA") 1 / 2 , X k = k — 1,2), $ 
becomes 

f°° 2tt 
$ = 27tto / exp (— toAcosIix) sinh^dx = -r- exp (— A) . (A15) 
Jo A 

From Eq. (A9)-(A10) and |(A15)|, it follows that 



rtf = 5gexp(-A) 



A 3 
2ttC 



{A 2 + 3A + 3}^^ + {A + l}^ 
A A 



exp(-A){A + l}A c 



(A16) 
(A17) 



On the other hand, the stress energy tensor of a perfect fluid is 

T a f 3 = (eo+Po)U a Uf 3 +p 5 a f 3 , (A18) 
while the numerical surface density of particles measured in the rest frame of the fluid S is related to via 

E = -<p«U a . (A19) 



Eq. (A 16) has clearly the form of Eq. (A18) when identifying 

U a = X a /X , 
2irC 

£q+Po = -^-exp(-A) {A 2 + 3A + 3} , 



Pa 



A 3 
2ttC 



exp (-A){A+ 1} 



(A20a) 
(A206) 

(A20c) 
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Eliminating A,C between Eqs. (Al7)-(A20c) yields the sought after equation of state 

E 



(A21) 



l+Po/E ' 

where po — E/A. This is the familiar gas law given that 1/A is the absolute temperature. The equation of 



state 



Eq. (A21) corresponds by construction to an isentropic flow. Indeed, using the conservation equations 



Eq. (Alia) dotted with U a , and Eq. (A116) together with U a U a = —1 yields after some algebra to the 
identity 



ds 



(A 2 + 3A + 3) 
A(A+1) 



1 

+ A 



dU a 
dx a 



1 dE 

E ds 



(A22) 



where d/ds = (dx a /dx a )(d/dx a ) is the covariant derivative following the stream lines. Writing the r.h.s of 
Eq. (A22) as an exact covariant derivative leads to the first integral 



d 
ds 

Given that po — S/A, it follows 



1 + A 



2 log (A) + log (1 + A)- log (E) 



0. 



( 1 + t) 



Po 

e exp vi + E/ Po ; E 



Const. 



(A23) 



(A24) 



which in the low temperature limit gives E 2 po = Const. This corresponds to the correct adiabatic index 
7 = (2 + 2)/2. 



B. Distribution functions for counter-rotating relativistic disks 

For counter rotating disks in the absence of gravomagnetic forces, the inversion for distribution functions 
given in section 4.1 can be carried in close analogy to the Newtonian method given by Pichon & Lynden-Bcll 
(1995) pQ ], and is sketched here. The relativistic flow is characterised by its stress energy tensor, T^ al3 \ 
which, in Vierbein components, reads 

From the geodesic equation it follows that 

W M = -1 =* P (fl) = ^[ £ 2 -fe 2 W7 # - ftt ] 1/2 , (B.2) 

■J la 



where the line element in the disk is taken to be 

ds 2 = - 7tt dt 2 + lRR dR 2 + lu d(f> 2 . (fl.3) 

Differentiating Eq. (B.2) implies 

dP {R) de dh . . 

— 77T- = 77TT , while dP w = , (BA) 

P {t) [e 2 -ft 2 W7#-7«f /2 



where the differentiation is done at constant R and h , and R and e respectively. Combinations of Eqs. (B.l) 
lead to 



A 



700 V2 ( T {tt) _ T {H) _ T (RR)) = f f f*(e,h)dhde ^ (fi ^ 



A = ,^ 2 TV* = I I f^ h ] dhd < . {BM) 

^/e A - h^tt/lH ~ Itt 
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Note that Eq. (B.ba) corresponds to Tolman's formula [p3| which yields be integration the total gravitating 



mass of the disk. Introducing 



K 2 = lu / ltt , e = i(e 2 -l), * = I(l_ 7tt ) an d £ ( e , h ) = LM1 j (R6) 



gives, for Eqs. (B.5): 



A=// Me,h)dhde and A = /" /" ^ Me,h)dhde 

[2(e + ^)-h 2 /K 2 ] 1/2 J J [2{e + iP)-h 2 /K 2 ] 1/2 



Eq. (g.7)| is formally identical to Eq.(2.9) of Pichon & Lynden-Bell [g0[, replacing (A, A, /*,7£, e, ^) by 
(i?S,i? 3 p0, f,R,e,ip). In fact (A, A, 7?., e, tends to (RT,, R 3 p^, f,R,e,ip) in the classical regime. 
Introducing the intermediate functions F(h,TZ) 1 which is chosen so that its moments satisfy Eqs. (-B.5), and 
following the classical prescription of Pichon & Lynden-Bell |^(| yields 

n r p (dP/dn) dn ., 7 (dP/dn) dn 

/(**) = = / ? T = ~~ / 7^ ^ 5 . 



7T 7 Jh 2 /n 2 - 2* - 2e 7T 7 Jh 2 /n 2 - 2* - 2e 
where lZ p (h, e), and lZ a (h, e) are respectively the apogee and perigee of the star with invariants (fe, e), and 



TZ e (h) is the inner radius of a star on a "parabolic" (zero e energy) orbit with momentum h. Eq. (B. 
provides direct and systematic means to construct families of distribution function for counter-rotating disks 
characterised by their stress-energy tensor in the plane. Note that in contrast the dragging of inertial frames 
requires to fix three moments of the velocity distribution to account for the given energy density, pressure 
law and rotation law. 

C. Distribution functions for rotating disks: Lorentzian in momentum 

A possible choice for F is a Lorenzian parametrised in width, p2, mean, p\ and amplitude, pq such that 

P (v ( f > ,R)dv < f > — PoP27r — ^E— y where p = an d v< f, = e" (x - /R c . {CI) 

P 2 2 + (P-Pi) Jl-v 2 



The explicit function of is found using dp = dv$/(l — v^fl 2 . In this instance, Eq. (4.12) becomes 

T<«> = ™*r -. * r = Pa , (C.2o) 



00 (]P2 2 + (P-Pl) 2 ) 

TW e- = ™»r / dp ^ = g (p , Pl , P2 ) , (C.2b) 

T (00) e 4„ = P0P2 f°° P^dp = Pa(p2 + P2 2 + Pi 2 ) ^ 2 ^ 

n J-oo {l+p 2 f (p 2 2 + (P ~ Plf) (1+P2) 2 +Pl 2 

where G(po,Pi,P2) is given by Eq. (C.4) bellow. The inversion of Eqs. (C.2) for (po, Pij P2), together with 
Eqs. (C.l), (4.14) and Eq. (4.3c) yields the distribution function of a Lorenzian flow compatible with the 
energy distribution imposed by the metric (is, £, uj) and the cut z — f(R). This inversion is always achievable 
provided the energy conditions Eq. (3.12) are satisfied. 

More generally, families of distributions functions corresponding to azimuthal momentum distribution 
such as 2p\pQ 7r _1 /(p 2 + (p — pi) 2 ) 2 , 8P2P0 (^ 7T )' 1 /{P2 + (J 3- Pi) 2 ) 3 ' ' ' are readily derived in a similar manner 
while differentiating the r.h.s. of Eq. (C.2) with respect to p\. These higher order distributions have all 
vanishing number of stars with azimuthal velocity close to the velocity of light. This parametrisation is 
illustrated for the Kerr disk on Fig. 7-4 
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C.l. The quadrature for G(po>PiiP2) 



Changing variables to u such that p = (u — l/u)/2 Eq. (C.2b) becomes 



Q{po,PuP2) = —— , 2 2 , , 2 — : — z 3 ( c - 3 ) 

7T Jo 4p 2 U 2 + (m 2 - 1 - 2p! u) 



This integral is then carried and yields 



Po {P2 +ipi) log ( -2ip 2 + 2pi - 2 J I + {-ip 2 +pif 



2 \fl + (~ip2 +Plf 7T 



Po (iP2 -Pi) log I -2ip 2 + 2pi + 2 Jl + (-ip 2 + Pif I Po(p2 ~ ipi) log ( 2ip 2 + 2pi - 2 v/l + (ip 2 + Pif 



2a/1 + Hp2+Pi) 2 ^ 2Jl + (ip 2 +Pi) 2 vr 



Po (ipi -Pa) log (2ip 2 + 2p x + 2 J 1 + (ip 2 + Pi f 



2\Jl + (ip2 +Pl) 2 7T 



(C.4) 



D. The relativistic epicyclic frequency 

The equation for the radial oscillations follows from the Lagrangian 



C = -de 2v [i -u(f>) 2 - e 2C -- 2v (1 + f' 2 ) R 2 - e" 2 " 2 i? 2 , (£U) 

where ( ) stands here for derivation with respect to the proper time r for the star describing its orbit. In 
Eq. (D.l), 4> and t are ignorable which leads to the invariants: 

dC , dC , . 

d<j) at 

The integral of the motion for the radial motion follows from Eqs. (D.2) and U^U^ = — 1, namely: 

R 2 



e 2(-2u ^ + ; ,2) ^2 + £_ {h _ euj) 2 _ e 2 e -2u = _j 



This equation, together with its total derivative with respect to proper time provides the radial equation of 
motion, having solved for the angular momentum h and the energy e of circular orbits as a function of radius 
when equating both R and R to zero, the relativistic generalisation of the classical epicyclic frequency is 
defined here to be the frequency of the oscillator calling back linear departure from circular orbits. Hence 
the equation for radial perturbation reads 

6R + k 2 SR = 0, (D.3) 

which gives for n 2 : 

n 2 =e- 2Z (Z" - 2Z' 2 ) - e 2 e~ 2L (L" - 2L' 2 ) - 

V ' V ' (DA) 

e 2K [ (h - 6lj) 2 (K" + 2K' 2 ) - e (ft - ew) {AK'J + u") + eV 2 ] , 
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where ( )' stands in this Appendix for dj dR = d/dR + f d/dz and L, K, and Z are given in terms of the 
potential £ and v by: 



L = C + logv / r + ?2, {D. ha) 

Z = C + log y/1 + /' 2 - I/, (D.56) 

K = 2v - ( - log v/1 + J' 2 - log R. (D.5c) 

Here e and /i are known functions of i? given by the roots of: 

R = e- 2i e 2 = (/i - e ^) 2 e 2Ar + e~ 2Z , (D.6a) 

R =0 =► e- 2Z Z' = e- 2i L'e 2 + (D.66) 



(/i-ew)e 2Ar [(/i-ew)7V'-ea;']. 



E. Correction to Lemos' Self-Similar Disk 

The result given in section 6.1 agrees with that of Lemos, having made the following agreed corrections to 
his solution: 

• Eq.(L3.14) should read: (k + n) cot ±(k + n)ir = —4irp„ r with a minus sign, because the integration of 
9 is carried downwards. This gives n + k > 1 for the disks with positive pressures. 

• Lemos gives the stress energy tensor integrated through the plane using coordinate increment rather 
than proper length. Our result gives the latter which has an extra factor exp(£ — v) and gives the 
physical energy tensor per unit area. 

• The exponent of cos[(n + k)ir/2] is 2n(k — n)/(n + k) 2 rather than 4n(l — 2n)/(k + n) 2 . 
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ABSTRACT 

Complete sequences of new analytic solutions of Einstein's equations which describe 
thin super massive disks are constructed. These solutions are derived geometrically. The 
identification of points across two symmetrical cuts through a vacuum solution of Einstein's 
equations defines the gradient discontinuity from which the properties of the disk can be 
deduced. The subset of possible cuts which lead to physical solutions is presented. At 
large distances, all these disks become Newtonian, but in their central regions they exhibit 
relativistic features such as velocities close that of light, and large redshifts. Sections with 
zero extrinsic curvature yield cold disks. Curved sections may induce disks which are stable 
against radial instability. The general counter rotating flat disk with planar pressure tensor 
is found. Owing to gravomagnetic forces, there is no systematic method of constructing 
vacuum stationary fields for which the non-diagonal component of the metric is a free 
parameter. However, all static vacuum solutions may be extended to fully stationary 
fields via simple algebraic transformations. Such disks can generate a great variety of 
different metrics including Kerr's metric with any ratio of a to m. A simple inversion 
formula is given which yields all distribution functions compatible with the characteristics 
of the flow, providing formally a complete description of the stellar dynamics of flattened 
relativistic disks. It is illustrated for the Kerr disk. 
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1. Introduction 

The general-relativistic theory of rapidly rotating objects is of great intrinsic interest and has potentially 
important applications in astrophysics. Following Einstein's early work on relativistic collisionless spherical 
shells, Fackerell (1968) [10], Ipser (1969) [15] and Ipser & Thorne (1968) [22] have developed the general 
theory of relativistic collisionless spherical equilibria. The stability of such bodies has been further expanded 
by Katz (1975) [11]. Here the complete dynamics of flat disk equilibria is developed and illustrated, extending 
these results to differentially rotating configurations with partial pressure support. A general inversion 
method for the corresponding distribution functions is presented, yielding a coherent model for the stellar 
dynamics of these disks. Differentially rotating flat disks in which centrifugal force almost balances gravity 
can also give rise to relatively long lived configurations with large binding energies. Such objects may 
therefore correspond to possible models for the latest stage of the collapse of a proto-quasar. Analytic 
calculations of their structure have been carried into the post Newtonian regime by Chandrasekhar. Strongly 
relativistic bodies have been studied numerically following the pioneering work on uniformly rotating cold 
disks by Bardeen & Wagoner [1]. An analytic vacuum solution to the Einstein equations, the Kerr metric [12], 
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which is asymptotically flat and has the general properties expected of an exterior metric of a rotating 
object, has been known for thirty years, but attempts to fit an interior solution to its exterior metric have 
been unsatisfactory. A method for generating families of self-gravitating rapidly rotating disks purely by 
geometrical methods is presented. Known vacuum solutions of Einstein's equations are used to produce the 
corresponding relativistic disks. One family includes interior solutions for the Kerr metric. 

1.1. Background 

In 1919 Weyl [26] and Levi-Civita gave a method for finding all solutions of the axially symmetric Einstein 
field equations after imposing the simplifying conditions that they describe a static vacuum field. Bicak, 
Lynden-Bell and Kalnajs (1993) [3]- here-after refered as BLK - derived complete solutions corresponding 
to counter rotating pressure-less axisymmetric disks by matching exterior solutions of the Weyl- Levi-Civita 
type on each side of the disk. For static metrics, the vacuum solution was obtained via line superposition 
of fictitious sources placed symmetrically on each side of the disk by analogy with the classical method of 
images. The corresponding solutions describe two counter-rotating stellar streams. This method will here be 
generalised, allowing true rotation and radial pressure support for solutions describing self-gravitating disks. 
These solutions provide physical sources (interior solutions) for stationary axisymmetric gravitational fields. 

1.1.1. Pressure Support via Curvature 

In Newtonian theory, once a solution for the vacuum field of Laplace's equation is known, it is straight- 
forward to build solutions to Poisson's equation for disks by using the method of mirror images. In general 
relativity, the overall picture is similar. In fact, it was pointed out by Morgan & Morgan (1969) [18] that for 
pressureless counter rotating disks in Weyl's co-ordinates, the Einstein equation reduces essentially to solving 
Laplace's equation. When pressure is present, no global set of Weyl co-ordinates exists. The condition for 
the existence of such a set of co-ordinates is that Trr + T zz = 0. On the disk itself, this is satisfied only 
if the radial eigenvalue of the pressure tensor, pr, vanishes (T zz vanishes provided the disk is thin). When 
Pr is non-zero, vacuum co-ordinates of Weyl's type still apply both above and below the disk, but as two 
separate co-ordinate patches. Consider global axisymmetric co-ordinates (Ft! ,<f>' , z') with z' = on the disk 
itself. For z' > 0, Weyl's co-ordinates R, <f>, z may be used. In terms of these variables, the upper boundary 
of the disk z' = + is mapped onto a given surface z = f(R). On the lower patch, by symmetry, the disk will 
be located on the surface z = —f(R). Thus in Weyl co-ordinates, the points z = f(R) on the upper surface 
of the disk must be identified with the points z = —f(R) on the lower surface. The intrinsic curvature of the 
two surfaces that are identified match by symmetry. The jump in the extrinsic curvature gives the surface 
distribution of stress-energy on the disk. A given metric of Weyl's form is a solution of the empty space 
Einstein equations, and does not give a complete specification of its sources (for example, a Schwarzschild 
metric of any given mass can be generated by a static spherical shell of any radius). With this method, 
all physical properties of the source are entirely characterised once it is specified that the source lie on the 
surface z = ±/(i?) and the corresponding extrinsic curvature (i.e. its relative layout within the embedding 
vacuum spacetime) is known. The freedom of choice of v in the Weyl metric corresponds to the freedom 
to choose the density of the disk as a function of radius. The supplementary degree of freedom involved in 
choosing the surface of section z = f(R) corresponds classically to the choice of the radial pressure profile. 
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1.1.2. Rotating Disks 



For rotating disks, the dragging of inertial frames induces strongly non-linear fields outside the disk which 
prohibit the construction of a vacuum solution by superposition of the line densities of fictitious sources, 
nevertheless, in practise quite a few vacuum stationary solutions have been given in the literature, the most 
famous being the Kerr metric. Hoenselaers, Kinnersley & Xanthopoulos (1979) [13,14](HKX hereafter) 
have also given a discrete method of generating rotating solutions from known static Weyl solutions. The 
disks are derived by taking a cut through such a field above all singularities or sources and identifying it 
with a symmetrical cut through a reflection of the source field. This method applies directly to the non- 
linear fields such as those generated by a rotating metric of the Weyl-Papapetrou type describing stationary 
axisymmetric vacuum solutions. The analogy with electromagnetism is then to consider the field associated 
with a known azimuthal vector potential A^, the analog of gt^/gtt- The jump in the tangential component 
of the corresponding 5-field gives the electric current in the plane, just as the jump in the t, 4> component 
of the extrinsic curvature gives the matter current. 

The procedure is then the following: for a given surface embedded in curved space-time and a given 
field outside the disk, the jump in the extrinsic curvature on each side of the mirror images of that surface 
is determined in order to specify the matter distribution within the disk. The formal relationship between 
extrinsic curvatures and surface energy tensor per unit area was introduced by Israel and is described in 
details by Misner Thorne Wheeler (1973) p 552. [19] Plane surfaces and non-rotating vacuum fields (i.e. the 
direct counterpart of the classical case) were investigated by BLK and led to pressureless counter-rotating 
disks. Any curved surface will therefore produce disks with some radial pressure support. Any vacuum 
metric with g t< f> ^ outside the disk will induce rotating disks. 



Z 




R 



Figure 1.1: Symbolic representation of a section z = ±/(i?) through the g t< f> field repre- 
senting the upper and lower surface of the disk embedded in two Weyl-Papapetrou fields. 
The isocontours of uj(R, z) correspond to a measure of the amount of rotation in the disk. 



4 



In section 2, the construction scheme for rotating disks with pressure support is presented. In section 3, 
the properties of these disks are analysed, and a reasonable cut, z = f(R), is suggested. Section 4 describes 
stellar equilibria which have the corresponding anisotropic stress energy tensors, while section 5 derives the 
properties of the limiting counter-rotating disks. This method is then first applied to warm counter rotating 
disks in section 6; the dominant energy condition given for Curzon disks by Chamorow el a/.(1987) [5], and 
Lemos ' solution for the self - similar Mestel disk with pressure support are recovered. In section 7, the Kerr 
disk is studied and the method for producing general HKX disks is sketched. 



2. Derivation 

The jump in extrinsic curvature of a given profile is calculated and related to the stress energy tensor 
of the matter distribution in the corresponding self gravitating relativistic disk. 

2.1. The Extrinsic Curvature 

The line element corresponding to an axisymmetric stationary vacuum gravitational field is given by 
the Weyl-Papapetrou (WP) metric 

d 8 2 =-e 2v {dt-ud<t>) 2 +e 2 <- 2v (^2 + rfz 2^ +i? 2 e -2, # 2 ; (21) 

where (R,<f>,z) are standard cylindrical co-ordinates, and v, ( and u> are functions of (R, z) onlygeo- 
metric units G = c = 1 are used throughout this paper. Note that this form is explicitly symmetric under 
simultaneous change of <f> and t. The vacuum field induces a natural metric on a 3-space-like z = f(R) 
surface 

dx a dx^ 

da 2 = g afj -^—^dx a dx l = g afj h^h^ dx a dx b , 

= jab dx a dx b , 

where {x a } are the co-ordinates on the embedded hypersurface z = f(R), {x a } = (t,R,<f>), and {x a } = 
(t, R, z, <f>) are Weyl's co-ordinates for the embedding spacetime. Here, j a i, stands for the embedded metric, 
and h" is given by 

"10 



dx° 
dx a 



10/' 
10 



(2.3) 



where hereafter a prime, ( )', represents the derivative with respect to R. The line element is then 

da 2 =-e 2v {dt-ud<t>) 2 + R 2 e- 2v d<t? + e 2 t- 2v (l + f 2 ) dR 2 . (2.4) 
Let the upward pointing normal to the surface z = f(R) be 

N fi = ^[z-f(R)] = (0,-fi,0,l). (2.5) 



The normalised normal is therefore 



N » - N » e^. (2.6) 



VN^W y/g"PNpN a y / TTr 2 ~ 
The extrinsic curvature tensor K is given in its covariant form by 

K ab =n»{d a K+T» ap h a a h(), (2.7) 
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where d a h^ = dh^/dx a . Eq. (2.6) and (2.3) together with the Christoffel symbols for the Weyl metric lead 
to the computation of K , the extrinsic curvature of that surface embedded in the Weyl Papapetrou metric: 



2C 



K 



RR 



V^Tr Vi + /' 



dN 



(i+n 



K, 



al>-(, 



K 



4>t 



4>4> 



R dN 



where the notation d/dN = (d/dz — f'd/dR) has been used. 
2.2. The Stress Energy Tensor 



d[v+±log(u)] 

dN 



dN 



(2.8a) 
(2.86) 
(2.8c) 
(2.8d) 



The stress energy tensor per unit surface t£ is the integral of of the stress energy tensor carried along 
the normal to the surface z = f(R). In a locally Minkowskian frame co-moving with the mean flow of the 
disk, the corresponding orthonormal tetrad is 



e (°W(i,o,-"), 

e^ 1 W"Vl + /' 2 (0,1,0), 
e^ = Re~ v (0,0,1), 



(2.9a) 
(2.96) 
(2.9c) 



so that 



d s 2 = V(a)(i) (e^dx 1 ) [e^dx 3 



with f](a)(b) = Diag(— 1, 1,1). In that frame, 



(2.10) 



-(«)(») 



£00 
Pn 

P<p 



(2.11) 



After a Lorentz transformation to a more general frame in which the flow is rotating with relative velocity 
V in the <j> direction, this stress becomes 



r(<0(») 



1 - V 2 



s+ P4> V 2 (p* + e)V 

(1 - V 2 ) p R 

(P4 + s)V P4> + zV 2 



(2.12) 
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2.3. The Discontinuity Equations 

Israel (1964) [16] has shown that Einstein's equations integrated through a given surface of discontinuity 
can be re-written in terms of the jump in extrinsic curvature, namely 



T a — — \K a — K a s a w = r a 

b o l Ix b Ix a °bl- — >->b ■ 
07T 



(2.13) 



where [ stands for ( ) taken on z = f(R) minus ( ) taken on z = —f(R). This tensor is known as 
the Lanczos tensor. 

In the tetrad frame (2.9), the Lanczos tensor given by Eqs. (2.13) and (2.8) reads 



£(«)(») 



4^1 + /' : 



l+f' 2 R dN v dN 2R 



dN 2R 



R 




l+f' 2 1 dN 



(2.14) 



Identifying the stress energy tensor with the tetrad Lanczos tensor £<»(') according to Eq. (2.13), 

and solving for p R ,p^, , e and V yields 



V 



Re 



-2v 



duj/dN 



f dv\ 

R +2 lm)- Q , 



47T y/1 + P 



Pr 



4^1 + /' : 



Q 
2 

Q 
2 



1 + /' 2 2R dN 
f" , /' , d(v-C) 



l + f' 2 2R 



dN 



-r_ 

4tt ^\ + f' 2 V R 



where 



Q 



R dN J R 2 \dN J 



(2.15a) 
(2.156) 
(2.15c) 
(2.15rf) 

(2.16) 



All quantities are to be evaluated along z = f(R). Eq. (2.15) gives the form of the most general solution to 
the relativistic rotating thin disk problem provided the expressions for e,p^,p R are physically acceptable. 



3. Physical properties of the warm disks 



Physical properties of interest for these disks are derived and related to the choice of profile z = f(R) 
compatible with their dynamical stability. 

Defining the circumferential radius R c , proper radial length R, and a synchronised proper time r* by 



R c = Rt~" , 

R= f + f' 2 eS- v dR, 



on 



e v [ 1 -w-^ ] dt 



(3.1a) 
(3.16) 

(3.1c) 



where the integral over dR is performed at z = f(R) and that over dt is performed along a given trajectory 
(T), the line element on the disk (2.4) reads 



2 i n2ji2 



la 2 = -ir 2 + dR + R\ 



(3.2) 
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For circular flows (^fr* = e"[l — u>d<f>/dt] dt), V may be re-written as 

Equation (3.3) is inverted to yield the angular velocity of the flow as measured at infinity: 

*= d ±= Ve " /R : . (3.4) 

dt l + ujVe v /R c v ; 

Similarly, the covariant specific angular momentum A of a given particle reads in terms of these variables 

h=^= 7t4>U t + 74>4> ^ = ^J==(R c V + u 1 e") , (3.5) 

where \i is the rest mass of the particle. The covariant specific energy e of that particle reads 

^-^ = - lt ^- lt y = - / J==. (3.6) 

which gives the central redshift, 1 + z c = exp {—v). The velocity of zero angular momentum observers (so 
called ZAMOs) follows from Eq. (3.5): 

cue" , x 

^=—5-- ( 3 - 7 ) 

K c 

With these observers, the cuts z = f(R) may be extended inside ergoregions where the dragging of inertial 
frames induces apparent superluminous motions as measured by locally static observers. The circumferential 
radius R z as measured by ZAMOs is Lorentz contracted with respect to R c , becoming 

R z = ^=R c [l-V^} 1 '\ (3.8) 
while the velocity flow measured by ZAMOs is 

V — V, 

The total angular momentum of the disk H follows from the asymptotic behaviour of the vacuum field u> 
at large radii, u> — ► —2H/r. Define the binding energy of the disk AE as the difference between the baryonic 
rest mass Mo and the total mass-energy of the disk M as measured from infinity given by v —> —M/r: 

AE = M - M. (3.10) 

M is also given by the Tolman formula, but corresponds, by construction, to the mass of the line source for 
the vacuum field. The baryonic rest mass can be expressed as 

M = j ^V^gU^dS^ E2wR e dR, (3.11) 

where the baryonic energy density £ can be related to the energy density e through some yet unspecified 
equation of state or the knowledge of a distribution function. 
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3.1. Constraints on the internal solutions 



The choice of the profile z = f(R) is open, provided that it leads to meaningful physical quantities. 
Indeed it must lead to solutions which satisfy £>0,p R >0,p^,>0, and!/ < 1 while the dominant energy 
condition implies also that e > p^ and e > p R . These translate into the following constraints on / 

e>0 => \jAl 2 - O 2 + (AT - 2Z) > ; (3.12a) 

0<V<1 => Af>0 and £>>0; (3.126) 
0< Pr p<e => 4Z Af > (O 2 + 4Z 2 ) andW>2Z; (3.12c) 

< P R < £ => /' < and - 2/7 R < \J M 2 - O 2 + (Af - 2Z) ; (3.12c?) 



given 



* s £ + 2 &. o s ^, »„ d z^fxry (3 . 13) 



R dN ' R dN ' \dN 1 + /■ 



The existence of a solution satisfying this set of constraints can be demonstrated as follows: in the limit 
of zero pressure and counter-rotation (i.e. O —> 0, Af —> 2dv/dz, and Z — ► d(/dz), any cut z = Const = 
b m satisfies Eqs. (3.12). By continuity, there exist solutions with proper rotation and partial pressure 
support. In practise, all solutions given in sections 6 and 7 satisfy the constraints (3.12). Note that in the 
limit of zero radial pressure, Eq. (3.126) implies d/dz [u> / R + exp 2v\ < 0. 



3.2. Ansatz for the profile 

In the following discussion, the section z = f(R) is chosen so that the corresponding radial pressure 
gradient balances a given fraction of the gravitational field which would have occurred had there been no 
radial pressure within the disk. This choice 'bootstraps' calculations for all relevant physical quantities in 
terms of a single degree of freedom (i.e. this fraction rj), rather than a complete functional. This gives 

Hk'-TTvr)^''^-'^ (3 - 14) 

On the r.h.s. of Eq. (3.14), p R is put to zero and V, p^,s are re-expressed in terms of (,i>, to via Eqs. (2.15) 
with /' = /" = and z = Const. On the l.h.s. of Eq. (3.14), p R is chosen according to Eq. (2.15c?). In 
practise, this Ansatz for / is a convenient way to investigate a parameter space which is likely to be stable 
with respect to ring formation as discussed in the next subsection. In principle, a cut, /, could be chosen so 
as to provide a closed bounded symmetrical surface containing all fictitious sources. Indeed, by symmetry 
no flux then crosses the z = plane beyond the cut, and the energy distribution is therefore bound to the 
edge of that surface: this corresponds to a finite pressure supported disk. It turns out that in general this 
choice is not compatible with all the constraints enumerated in the previous section. More specifically, the 
positivity of p^ fails for all finite disk models constructed. The Ansatz described above gives for instance an 
upper bound on the height of that surface when assuming that all the support is provided by radial pressure. 
This height is in turn not compatible with positive azimuthal pressure at all radii. 
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3.3. Stability 

To what extent are the equilibria studied in the previous sections likely to be stable under the action of 
disturbances? The basic instabilities can be categorised as follows: 

dynamical instabilities, which are intrinsic to the dynamical parameters of the disks, grow on an orbital 
time scale, and typically have drastic effects on the structure of the system. 

secular instabilities, which arise owing to dissipative mechanisms such as viscosity or gravitational 
radiation, grow on a time scale which depends on the strength of the dissipative mechanism involved, and 
slowly drive the system along a sequence of dynamical equilibria. 

Amongst dynamical instabilities, kinematical instabilities correspond to the instability of circular orbits 
to small perturbations, and collective instabilities arises from the formation of growing waves triggered by 
the self-gravity of the disk. Rings, for instance, will be generated spontaneously in the disk if the local radial 
pressure is insufficient to counteract the self-gravity of small density enhancements. Even dynamically stable 
non-axisymmetric modes may drive the system away from its equilibrium by radiation of gravitational waves 
which will slowly remove angular momentum from the disk. For gaseous disks, viscosity and photon pressure 
will affect the equilibria. Radiative emission may disrupt or broaden the disk if the radiation pressure 
exceeds the Eddington limit. The energy loss by viscosity may induce a radial flow in the disk. However, 
the disks discussed in this paper have anisotropic pressures inappropriate for gaseous disks and the accurate 
description of the latter two processes requires some prescription for the dissipative processes in the gas. 
The scope of this section is therefore restricted to a simple analysis of the dynamical instabilities. 

Turning briefly to the corresponding Newtonian problem, Toomre (1963) [24] gave the local criterion 
for radial collective instability of stellar disks, 

a R > , (3.15) 

K 

where So is the local surface density, u 2 R the radial velocity dispersion, and k the epicyclic frequency of the 
unperturbed stars. This criterion is derived from the first critical growing mode of the dispersion relation 
for radial waves. The corresponding critical wavelength is Jeans length ~ 2<j 2 r / GY*®. 

For the relativistic disks described in this paper, spacetime is locally flat, which suggests a direct 
translation of Eq. (3.15) term by term. The constraints that stability against ring formation places on these 
models will then be addressed at least qualitatively via the Newtonian approach. The proper relativistic 
analysis is left to further investigation. The relativistic surface density generalising So is taken to be the 
co-moving energy density s given by Eq. (2.156). The radial velocity dispersion a 2 R is approximated by p R /s. 
The epicyclic frequency is calculated in the appendix. Putting Eqs. (DA) and (-D.5) into Eq. (3.15) give 
another constraint on / for local radial stability of these disks. Note that the kinematical stability of circular 
orbits follows from Eq. (D.3) by requiring k 2 to be positive. 



4. Stellar dynamical equilibria for rotating super-massive disks 

The method described in section 2 will generally induce rotating disks with anisotropic pressures (p^ ^ 
Pr). These objects may therefore be described in the context of stellar dynamics. It should be then checked 
that there exist a stellar equilibria compatible with a given vacuum field and a given cut z = f(R). 
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4.1. Distribution functions for rotating super-massive disks 

Here a general procedure to derive all distribution functions corresponding to a specified stress energy tensor 
is presented for disks with non-zero mean rotation (a more direct derivation for counter-rotating disk is given 
in Appendix B). The detailed description of the dynamics of the disk follows. 
In Vierbein components, the stress energy tensor reads 

Ti"P) = J J fAc,h)P^PW dp( R) dpW ^ (4 1} 

where /*(e, h) is the distribution of stars at position R, <j> with momentum P^ R \ P^\ For a stationary disk, 
it is a function of the invariant of the motion e, h. Now for the line element (2.4): 

da 2 = - e 2u (dt - uj d<j>f + t ic --' lv (l + /' 2 ) dR 2 + R 2 t~ 2v dcf, 2 , (4.2) 

the Vierbein momenta read 

P W = e v R^ih-Loe) , (4.3a) 
=e~ v e, (4.36) 

p(R) = [ £ 2 e -2, _ e 2u R -2 {h _ w f) 2 _ A 1/2 _ (4 3c) 



Calling x = h/e and $ = 1/e 2 , Eqs. (4.3) becomes 

pW =e v R- 1 ^- 1 ! 2 (x-lo) , (4.4a) 

p(t) =e - v 0-1/2 ^ ( 4 4& ) 

= i)- 1 ' 2 \e~ 2v - e 2v R- 2 (x - uj f - tfl ^ . (4.4c) 



In terms of these new variables, the integral element dP^dP^ b 

dp(4>) Qp(R) 



d pW d p( R ) 



dx dd 



d X d*= . o-"*- 1 *- 2 **** (4.5) 

2\ e~ 2v - e 2v R- 2 ( X - u>) 2 - 



Given Eq. (4.5) and Eq. (4.36), Eq. (4.1) may be rewritten as 

T («P) = f f„(a)„(0) k(^h)R-^-*l 2 dxd{) 



2\le- 2v - e 2v R- 2 ( X - uj) 2 - d 



In particular, 



RJ <u) e 2 V = f f U(e,h)#- 5 / 2 d X d# 

' ' 2^e~ 2v -e 2v R- 2 ( X - oj) 2 - d 
R 2 Tm = f / (x _ w) Ui^^dxM 



2Je~ 2 » - t 2 " R- 2 (x - to) 2 - d 



R*T^e- 2 » = [ [( X -.f , f*^)*-" 2 '*'* 

J J n - — '>,, '>,, n—'> . .\2 



2^e~ 2v - t 2 " R- 2 (x -u) -d 
Note that given Eqs. (4.7), T^ RR ^ follows from the equation of radial support. Note also that 



(4.6) 



(4.7a) 
(4.76) 
(4.7c) 



' ' 2\j e~ 2v - e 2v R~ 2 ( X - oj) 2 - d ' 
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or equivalently, 

T (U) _ T (H) _ rp(RK) = [ [ f*( € > h ) dp W d p(K) 



(4.9) 



Equation (4.9) yields the total gravitating mass. It is known as Tolman's formulae and can be derived 

directly from the relation P^P^ = —1 and Eq. (4.1). 

Let 

F(x,R) = \ I , f^- 5,2M J f-ixW* , (4 . 10) 

2J ^ e -2,_ e 2, R -2 ( X _ w ) 2 -^ { 2v^? 

where /* (x,$) = U (e,/i)^" 5/2 , and y (R, X ) = e~ 2v - e 2v RT 2 ( X - uj f . (4.11) 
The set of Eqs. (4.7) becomes 



RT (t ^e 2v = J F( X ,R) d X = Re~ 2v J F (^, R) dv^ (4.12a) 

R 2 T (t ^> = j F(x,R){x-u)d X = R 2 e~ Av j v <p F(v <p ,R)dv <p , (4.126) 
R 3 T (H) e -2, = l F {x,R) {x-uf d X =R 3 e~ 6v IvlFiy^R) dv^ . (4.12c) 



The r.A.sofEqs. (4.12) was re-expressed convinienlty in terms of the new variable = e" (x~ u>)/R c , and 
the new function F(v^, R) = F(x = R c v^e~ v +ui, R) in order to stress the analogy with the classical identities. 
Note that in the limit of zero pressures, reduces to V , the geodesic velocity of stars on circular orbits. 
In Eq. (4.12), the dragging of inertial frames requires to fix simultaneously three moments of the velocity 
distribution to account for the given energy density, pressures and rotation law. This third constraint does 
not hold in the corresponding Newtonian problem, nor for static relativistic disks as described in appendix B. 
Any positive function F(x, R) satisfying the moment constraints Eqs. (4.12) corresponds to a possible choice. 
For instance the parametrisation is carried in Appendix C for F distributions corresponding to powers of 
Lorentzians in azimuthal momenta. Let F(x,y) = F(x,R), where y(R,x) is given by Eq. (4.11). Written 
in terms of F , the integral equation Eq. (4.10) is solved for /* by an Abel transform 

1 X 

The latter integration may be carried in R space and yields 

* / u\ 2 f f dF \ dR 2 f f dF \ dR 

where p( R ) is given by Eq. (4.3c) as a function of R, e, and h. Here F is chosen to satisfy Eq. (4.12) which 
specifies the stress energy components T^"^ . The integration limits correspond to two branches: the first 
branch is to be carried between lZ e (h), the inner radius of a star on a "parabolic" orbit with momentum 
h, and lZ p (h,e), the perigee of a star with invariants (h,e); the second branch contributes negatively to 
Eq. (4.14) and corresponds to radii larger than 1Z a {h,e), the apogee of a star with invariants (h,e). Note 
that the derivative in Eq. (4.14) is taken at constant reduced momentum x = h/e. This inversion procedure 
is illustrated on Fig. 7.5. All properties of the flow follow in turn from /*(e, h). For instance, the rest mass 
surface density, S, may be evaluated from the detailed "microscopic" behaviour of the stars, leading to an 
estimate of the binding energy of the stellar cluster. 
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R 



Figure 4.1: sketch for the effective potential Y(R) as a function of radius. The horizontal 
line corresponds to the energy level 1/e 2 , the absisse line corresponds to the escape energy 
1/e 2 = 1. The ordinate axis corresponds to lZ e (h,e), the inner radius of a star on a 
"parabolic" orbit with momentum h, the left vertical line corresponds to lZ p (h,e), the 
perigee of a star with invariants (ft,e); the right vertical line corresponds to 1Z a {h,e), 
the apogee of a star with invariants (/i,e). The area between the curve Y = Y(R) and 
the lines Y = 1/e 2 and R = 1Z e are shaded, defining three regions from left to right. 
The middle region does not contribute to Eq. (4.14). The equation Y = 1/e 2 has two 
roots corresponding to lZ p (h, e) and 1Z a {h, e), while Y = 1 has two roots corresponding to 
infinity and lZ e . The sign of the slope of Y(R) gives the sign of the contribution for each 
branch in Eq. (4.14). 

4.2. A simple equation of state 

An equation of state for these rotating disks with planar anisotropic pressure tensors is alternatively found 
directly while assuming (somewhat arbitrarily) that the fluid corresponds to the superposition of two isotropic 
flows going in opposite directions. In other words, the anisotropy of the pressures measured in the frame 
co-moving with the mean flow V is itself induced by the counter rotation of two isotropic streams. (These 
counter rotating streams are described in more details in the next section.) In each stream, it is assumed 
that the pressure is isotropic and that the energy is exchanged adiabatically between each volume element. 
The detailed derivation of the corresponding relationship between the surface density and the pressure of 
each stream is given in Appendix A. The above set of assumption yields the following prescription for S : 

^ £ -Pf -P R + v^g -P<t> -Ph)(£ -P<t> + 3p R ) 

where Vo = \/ (p$ — P R )/{£ + P R ) is the counter rotating velocity measured in the frame co-moving with V 
which induces ^ p R . In the classical limit, E — ► (s — P4>)/\/l — V 2 . The binding energy of these rotating 
disks is computed from Eq. (4.15) together with Eqs. (3.11) and (2.15). 
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5. Application: warm counter-rotating disks 



For simplicity, warm counter-rotating solutions are presented first, avoiding the non-linearities induced 
by the dragging of inertial frames. These solutions generalise those of BLK, while implementing partial 
pressure support within the disk. Formally this is achieved by putting uj(R, z) identically to zero in Eq. (4.2); 
the metric for the axisymmetric static vacuum solutions given by Weyl is then recovered: 



where the functions £(R, z) and v(R, z) are generally of the form (cf. Eq. (BLK 2.21)) 

W (6) db 



R 2 + (\z\ + bf 
W(b 1 )W(b 2 )Z(b 1 ,b 2 ) db l db 2 , 



with Z given by 



Z 



2rir 2 



1 



r 2 - ri 
b 2 - bi 



> 0. 



(5.1) 

(5.2a) 
(5.26) 

(5.3) 



Here b\ and b 2 are the distances of two points below the disk's center and r\ and r 2 are the distances 
measured from these points to a point above the disk. Below the disk, Z is given by reflection with respect 
to the plane z = 0. When the line density of fictitious sources is of the form W(b) oc b~ m (Kuzmin-Toomre), 
or W(b) oc (*>( m ) (6) (Kalnajs-Mestel), the corresponding functions (,i> have explicitly been given in Bicak, 
Lynden-Bell & Pichon (1994) [4]. For the metric of Eq. (5.1), the Lanczos tensor (2.14) reads 



£(*)(*) 



47ryiT7 



72 



_X_ £ d(2 V -C) 
1 + f' 2 T R T dN 



-JL 



1+/' 



aw 



(5.4) 



Consider again counter-rotating disks made of two equal streams of stars circulating in opposite direc- 
tions around the disk center. The stress energy tensor is then the sum of the stress energy tensor of each 
stream: 



r («)(») 



e 





" 


2 




_0 


Pr 






P<P _ 


~ l-Vg 



s + Po Vg 

(1 - ^ 2 ) Po 

Po - 








(5.5) 



the pressure po in each stream being chosen to be isotropic in the plane. Subscript ( )o represents quantities 
measured for one stream. Expressions for s,p^, and p R follow from identifying Eqs. (5.5) and (5.4) according 
to Eq. (2.13). Solving for So, Vo, andpo in Eq. (5.5), given Eq. (5.4) and (2.13) gives 



Po 



47T y/1 + P 



f" 



i + /' 2 
-/'" 



dN 



87T^1 + f' 2 [ R . 

0L _ f" 

dN 1 + f' 2 + 



1 + 



72 



dN 



(5.6a) 
(5.66) 
(5.6c) 



The above solution provides the most general counter rotating disk model with pressure support. Indeed, 
any physical static disk will be characterised entirely by its vacuum field v and its radial pressure profile 
p R , which in turn defines W(b) and f(R) uniquely according to Eqs. (5.66) and (BLK 2.33). The other 
properties of the disk are then readily derived. 
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The angular frequency and angular momentum of these disks follow from Eq. (3.5), (3.4) on putting u> 
to zero and V to Vq. The epicyclic frequency k given by Eq. (DA) may be recast as 



2 _ t v dh 2 

which relates closely to the classical expression k 2 = R~ 3 dh 2 / dR. 

The Ansatz given by Eq. (3.14) for z = f(R) becomes, after the substitutions e —> So, V 



dv dv 



dv_ 1 5C 

dz 2 dz 
dv 



R 



dR 



1 



(5.7) 

V , p<p -> 0, 
(5.8a) 

(5.86) 



where z is to be evaluated at 6. Equation (5.86) follows from Eq. (5.8a) given the (^) component of the 
Einstein equation outside the disk. 

The equation of state for a relativistic isentropic 2D flow of counter-rotating identical particles is derived 
in appendix A while relating the perfect fluid stress energy tensor of the flow to the most probable distribution 
function which maximises the Boltzmann entropy. It reads 



s = 2p + 



Solving for £ gives 



Po 



(5.9) 



(5.10) 



which in the classical limit gives £ — ► Eq — po- The binding energy of these counter rotating disks is computed 
here from Eq. (5.10) together with Eqs. (3.11) and (5.6). Alternatively, it could be computed by constructing 
distribution functions using the inversion described in Appendix B. 



6. Examples of warm counter-rotating disks 



6.1. Example 1 : The Self-similar disk 



The symmetry of the self-similar disk allows one to reduce the partial differential equations corresponding 
to Einstein's field equations to an ordinary differential equation with respect to the only free parameter 
9 = arctan(i?/z) (Lemos, 89 [17]). Lemos' solution may be recovered and corrected by the method presented 
in sections 2 and 4. Weyl's metric in spherical polar co-ordinates (p, x,<i>) defined in terms of (R,z,<f>) by 
p = ^/(R 2 + z 2 ) and cot(x) = z/R is 



da 2 = -e 2v dr 2 + p 2 \e~ 2v sin 2 (*) d<f> 2 + e 2 <" 2 " ( ^ + d X 

L \p 

A self-similarity argument led Lemos to the metric 
ds 2 



-r 2n e N dt 2 



2P — N 7,2 . Z-N 

e d<p + e , ^ 2 



dr' 



where P , N , and Z are given by 



P (6) = log [sin (k + n) 6], 
An 

N (6) = log [cos (A; + n) 6/2], 



Z(0) 



k + n 
8n 2 
(k + nj' 



■ log [cos (k + n) 0/2] + 2 log (k + n) 



(6.1) 

(6.2) 

(6.3a) 
(6.36) 

(6.3c) 
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Here the <^'s in both metrics have already been identified since they should in both case vary uniformly in 
the range [0, 2ir[. The rest of the identification between the two metrics follows from the transformation 



x = (n + k) 
P = r n+k , 

T = t, 



with the result that v and £ take the form: 



2v = N + 2nlog(r) , 
2( = Z - 2 log (n + k) . 

This can be written in terms of the new variables given in Eqs. (6.4) as 

An 2 



c 



(k + n y 

2n 
(k + n) 2 



■log (f 
■ loe 



■ log (p) 



(6.4a) 
(6.46) 
(6.4c) 



(6.5a) 
(6.56) 



(6.6a) 
(6.66) 



The procedure described earlier may now be applied. In order to preserve the self-similarity of the solution 
and match Lemos' solution, f(R) is chosen so that 



/' (R) = const. = cot (rj) say. 



Then 



_5_ = fd_ _ ,d_ 

ON ~ \dz ■'dR 



1 



d 



p sin (rj) drj 



(6.7) 
(6.8) 



Putting Eqs. (6.6) and Eq. (6.8) into Eq. (2.15) gives 



e 

Pr 

,P4> 



2n(k — n) 

r ~ \ mi ( & +«) 2 

cos — 

4tt L V2, 



cot (rj) 



2n k 



tan (?j/2) 



2n^ 



(k+n)- 



(k+n) 
- cot (fj) 

■tan (77/2) 



(6.9) 



with rj to be evaluated at (n + k) tt/2. The above equations are in accordance with the solution found by 
Lemos by directly integrating the Einstein field equations. 



6.2. Example 2: The Curzon disk 



The Curzon disk is characterised by the following pair of Weyl functions: 

v = -a/r, ( = -a 2 R 2 / (2 r 4 ) , 



(6.10) 



with r = ^JR 2 + f 2 , given the dimensionless parameter a = M/b (recall that G = c = 1), all lengths being 
expressed in units of 6. This disk is the simplest example of the Kuzmin-Toomre family, where W(b) oc 6(b). 
It also corresponds to the building block of the expansion given in Eqs. (5.2). Equations (5.4) and (5.5) then 
imply 



e 

Pr 

,P4> 



4tt6^1 + f' 2 



+ £ + £ {2/ (r 3 - a R 2 ) + [a (R 2 - f 2 ) - 2r 3 ] Rf) 



1+/' 



T I^ + fi[2R 2 f-(R 2 -f 2 )Rf] 



(6.11) 




Figure 6.1: the "Toomre Q number" Q = <t r k/(3.36£o) is plotted here as a function of 
circumferential radius R c for different compactness parameter a = b/M when rj = 1/25. 

The weak energy conditions e > 0, e + p R > 0, and £ + > which follow are in agreement with 
those found by Chamorow et a./. (1987) [5]. Their solution was derived by direct integration of Einstein's 
static equations using the harmonic properties of the supplementary unknown needed to avoid the patch of 
Weyl metric above and below the disk. The cut z = f(R) corresponds to the imaginary part of the complex 
analytic function, the existence of which follows from the harmonicity of that new function. 

The Ansatz (3.14) leads here to the cut 



/(B) 



bm 2 2bm 3 (Ab 2 + 7R 2 ) 

4(b 2 + R 2 ) ~ !05(& 2 + i? 2 ) 5/2 



(6.12) 



In Fig. 6.1, the ratio of the binding energy of these Curzon disks (which is derived from Eqs. (3.11), 
(5. 10), (6. 11) and (6.12)) over the corresponding rest mass Mo is plotted with respect to the compactness 
parameter b/Mo for different ratios po/so measured at the origin. The relative binding energy decreases in 
the most compact configurations because these contain too many unstable orbits (k 2 < 0). A maximum 
ratio of about 5 % is reached. 
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Figure 6.2: the binding energy over the rest mass AE/Mo for Curzon disks, plotted 
against the compactness parameter a = b/M for different ratios po/so fixed by rj = 
0.01,0.16, 0.7. 

7. Examples of rotating disks and an internal solution for the Kerr metric 

The problem of constructing disks with proper rotation and partial pressure support is more complicated 
but physically more appealing than that of counter rotating solutions. It is now illustrated on an internal 
solution for the Kerr metric, and a generalisation of the solutions constructed by BLK and BLP to warm 
rotating solutions is sketched. 

7.1. A Kerr internal solution 

The functions (v, (,uj) for the Ker r metric in Weyl-Papapetrou form expressed in terms of spheroidal co- 
ordinates (R = s^(x 2 =p 1)(1 — t/ 2 ), z = sxy) are given by 



,2k 



-m 2 + s 2 x 2 + a 2 y 2 

2 ,.2 



,2C 



and 



where s 



(m + s x) + a 2 y 

—m 2 + s 2 x 2 + a 2 y 2 
s 2 (x 2 =p y 2 ) 

2am (m + s x) (l — t/ 2 ) 
— m 2 + s 2 x 2 + a 2 y 2 



(7.1) 
(7.2) 

(7.3) 



(±(m 2 



J )) 1/2 . H ere ± corresponds to the choice of prolate (+), or oblate (— ) spheroidal 



co-ordinates corresponding to the cases a < m and a > m respectively. In this set of co-ordinates, the 
prescription given in section 2 leads to completely algebraic solutions. The normal derivative to the surface 
z = f(R) reads (when a > m) 



(8/dN) = N (x, y) (d/dx) + N (y, x) (d/dy) , 



(7.4) 
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where 

N(x,y) = (x 2 -l) [y/s + x (y 2 - l) f'/R] /(x 2 - y 2 ) . 

Differentiating Eqs. (7.1) through (7.3) with respect to x and y together with Eq. (7.4) leads via Eqs. (2.15) to 
all physical characteristics of the Kerr disk in terms of the Kerr metric parameters m and a, and the function 
z = f(R) which must be chosen so as to provide relevant pressures and energy distributions according to 
Eqs. (3.12). (Bicak & Ledvinka (1993) [2] have constructed independently a cold Kerr solution when a < m.) 
On the axis R = 0, / = 6, and /" = — c, while Eqs. (3.12) imply /' = at (x = b/s, y = 1). This in turn 
implies 




£ / _J_ — fjl 2 

(P4>)o = (Pr)o = ,., , , t , ___,2 > ( 7 - 5a ) 



(6 + m) 




b 2 



2 



(6 + m) 2 



2m I (6 + m) 2 



(a 2 + (6 + m) 2 ) (a 2 + b 2 - m 2 ) 



(7.56) 



where ( )o stands for ( ) taken at R = 0. The constraint that all physical quantities should remain positive 
implies 



b > \J m 2 — a 2 , (7.6a) 

(b + m) 2 — a 2 , 
c<2m- v J . (7.66) 



(a 2 + (6 + m) 2 ) (a 2 + 6 2 - m 2 ) 



Equation (7.6a) is the obvious requirement that the surface of section should not enter the horizon of the 
fictitious source. Similarly, ergoregions will arise when the cut z = f(R) enter the torus 

z\ 2 ^/i? a ^ lm\ 2 /lm\ 2 (7 7) 

ml \m m 2 a J \2 a J 

The analysis may be extended beyond the ergoregion via the ZAMO frame. The characteristics of these 
frames are given by Eq. (3.7), (3.8) and (3.9). The central redshift 

l + z e = e -= ^ + y + i , (7.8) 

b 2 + a 2 — m 2 



can be made very large for such models constructed when 6 — ► \/m 2 — a 2 . 

On Fig. 7.1, Fig. 7.2 and Fig. 7.3, some characteristics of these Kerr disks with both a > m and 
a < m are illustrated. For simplicity, pressure is implemented via the cut (6.12). These disks present 
anisotropy of the planar pressure tensor (p^ ^ p R ). The anisotropy follows from the properties of the u> 
vacuum field which gives rise to the circular velocity curve. Indeed, in the outer part of the disk, the specific 
angular momentum h tends asymptotically to a constant. The corresponding centrifugal force is therefore 
insufficient to counter-balance the field generated by the v function. As this construction scheme generates 
such equilibria, the system compensates by increasing its azimuthal pressure away from isotropy. When 
a > m, and rj = the height of the critical cut corresponding to the last disk with positive pressures 
everywhere scales like a (a 2 — m 2 ft' 2 in the range 2 < a/m < 20. Note that this limit is above that of the 
highly relativistic motions which only occurs when 6 tends asymptotically to the Kerr horizon (a 2 — m 2 ) 1 / 2 . 
The binding energy which follows Eq. (4.15) reaches values as high as one tenth of the rest mass for the most 
compact configurations. The inversion method described in section 4.1 was carried for the Kerr disks and 
yields distribution functions which characterise completely the equilibria as illustrated on Fig. 7.5. It was 
assumed that the velocity distribution of stars at radius R was a modified squared Lorentzian, as discussed 
in Appendix C and shown on Fig. 7.4. Note that disks with low rotation support present in their outer 
parts two distinct maxima, one of which is counter rotating, in agreement with the assumptions yielding the 
equation of state Eq. (4.15). 
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Figure 7.1: the azimuthal velocity (a), the radial pressure (b), the ratio of radial over 
azimuthal pressure and the relative fraction of Kinetic Support (d) namely sV 2 /(sV 2 + 
P<P +Pr) f° r the a/m = 0.5 Kerr disk as a function of circumferential radius R c , for a class 
of solutions with decreasing potential compactness, b/m = 1.1, 1.2 ••• 1.7, and relative 
radial pressure support rj/m = 1/5 + 2(b/m — l)/5 . 

7.2. Other rotating solutions 



The rotating disk models given in section 2 require prior knowledge of the corresponding vacuum solu- 
tions. However, while studying the symmetry group of the stationary axially symmetric Einstein Maxwell 
equations, Hoenselaers, Kinnersley & Xanthopoulos (1979) [13,14] found a method of generating systemat- 
ically complete families of stationary Weyl-Papapetrou vacuum metrics from known static Weyl solutions. 
Besides, the decomposition of the Weyl potential vq into line densities provides a direct and compact method 
of finding solutions to the static field. The N+l rank zero HKX transform is defined as follows: let vq and 
Co be the seed Weyl functions given by Eqs. (5.2). HKX define the N x N matrix T parameterised by N 
twist parameters and N poles, a^, k = I, ■ ■ - N as 



'k,k' 



rk 



rk ~ r k ' 



± 1 



dk — O-k 1 



(7.9) 



with = r(aj.) = yj ' R 2 + (z — aj,) 2 and = j3(a]~) , j3 being a function satisfying (7.11a) below. Defining 
r = Diag(rj,), they introduce the auxiliary functions 

D ± = |l + r ± | , and L± = 2D ± tr [(l + T ± ) _1 T ± r] , 

w± = D+ ± e 2vo D- , and M± = zuu± + 2(L + T e 2v °L_) , 
which lead to the Weyl-Papapetrou potentials 



~D~' 

D+ 



(7.10a) 
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Figure 7.2: as in Fig. 7.1 with a/m = 0.95, rj/m = 1/15 + (6 — l)/5. Note the large 
relative kinetic support. 



e 2 < =ke 2 <°sR[D-D+* 
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l w +| 2 - l w -| 2 



(7.106) 
(7.10c) 



where 3£ and 9 stand respectively for the real part and the imaginary part of the argument and ( )* 
represents the complex conjugate of ( ). A; is a constant of integration which is fixed by the boundary 
conditions at infinity. The two yet unspecified function [3 and w satisfy: 



V/3 (a k ) = — ([z - a k ] V + Re^, x V) v , 
rk 



(7.11a) 
(7.116) 



where V = (d/dR, d/dz) and is the unit vector in the <f> direction. The linearity of these equations 
suggests again that solutions should be sought in terms of the line density, W(b), characterising vq, namely 



P(a k ) 



W{b) 



R 2 + (z-b-a k y 
R 2 + (z- b f 



db, 



W(b)(b-z) 
R 2 + (z- b f 



db. 



(7.12a) 
(7.126) 



It is therefore a matter of algebraic substitution to apply the above procedure and construct all non-linear 
stationary vacuum fields of the Papapetrou form from static Weyl fields. Following the prescription described 
in section 2, the corresponding disk in real rotation is then constructed. The requirement for the physical 
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Figure 7.3: the azimuthal velocity (a), the azimuthal pressure (b), the relative fraction 
of Kinetic Support (c), namely sV 2 /(sV 2 +£></>) and the angular momentum (d) for the 
a/m = 10 zero radial pressure Kerr disk as a function of circumferential radius R c , for a 
class of solutions with decreasing potential compactness b/m = 210, 260 • • -510. The first 
curve corresponds to a cut which induces negative azimuthal pressure (tensions). 




Figure 7.4: Parametrization of the number of star with velocity at radius R as a 
squared Lorentzian (as discussed in Appendix C) for the Kerr disk. The corresponding 
radius is appended on each curve. The left panel correspond to a = 0.2 while the right 
panel corresponds to a = 0.8. The pressure law is given in the caption of Fig. 7.5 where 
the corresponding distribution function is illustrated. 



source to be a disk is in effect a less stringent condition on the regularity of the vacuum metric in the 
neighbourhood of its singularity because only the half space which does not contain these singularities is 
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Figure 7.5: the isocontours of the distribution function for the Kerr disk as a function 
of the reduced momentum x = h/e, and the relative eccentricity of the orbit e c = (e — 
e x )/(l — e x ) where e x is the energy of circular orbits with reduced momentum x- The 
four panels correspond to a = 0,0.1,0.2,0.3,0.4 from left to right and top to bottom, 
while the pressure law was chosen to correspond to the cut f(R) = — 4(4 + 7i? 2 )/(525(l + 
i? 2 ) 5 / 2 ) + (16 + 15i? 2 )/(10(l + R 2 )) corresponding to a tepid disk. The figures are centered 
on the inner parts of the disk to illustrate the relative shift of the maximum number of 
stars towards larger angular momentum for faster rotating disks. The inversion assumed 
a squared Lorentzian for the distribution in angular velocity as discussed in Appendix C 
and illustrated on Fig. 7.4. 



physically meaningful. A suitable Ehlers transformation on D-/D+ ensures asymptotic flatness at large 
distance from the source. 

To illustrate this prescription consider the vacuum field given by Yamazaki (1981) [27]. This field arises 
from 2 HKX rank zero transformation on the Zipoy-Voorhess [25] static metric given by 



Co =^ 2 log 



(x + iy 

(x 2 - y 2 ) ' 



(7.13a) 
(7.136) 
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Figure 7.6: the azimuthal velocity for the aj Ky zero radial pressure Yamazaki disk as 
a function of circumferential radius R c , for a class of solutions with decreasing potential 
compactness bjny = 3.6, 3.8, • • • 5. Ky is the natural unit of distance in these disks as 
defined by Yamazaki's Eq. (Y-2) [27]. 



in prolate spheroidal co-ordinates. It corresponds to a uniform line density between ±6. The functions [3 
and w follow from Eqs. (7.12). Choosing poles a^, k = 1, 2 at the end of the rod ±6, leads to 



/?(a±) 



±<51oa 



26y. 



(7.14a) 
(7.146) 



Equations (7.10) together with (7.13) and (7.11) characterize completely the three Weyl-Papapetrou metric 
functions v, (, andw; all physical properties of the corresponding disks follow. For instance, Fig. 7.6 gives the 
zero radial pressure velocity curve of the Yamazaki disk spanning from the Schwarzschild (8 = 1/2) metric 
when u\ = «2- 

More generally, the extension of this work to the construction of disks with planar isotropic pressures 
should be possible by requiring that Eq. (2.15c) is identically equal to Eq. (2.15c?) given Eq. (7.10). 



8. Conclusion 



The general counter rotating disk with partial pressure support has been presented. The suggested 
method for implementing radial pressure support also applies to the construction of stationary axisymmetric 
disk with rotation but requires prior knowledge of the corresponding vacuum solution. In this manner, a 
disk-like source for the Kerr field has been constructed orbit by orbit. The corresponding distribution reduces 
to a Keplerian flow in the outer part of the disk and presents strong relativistic features in the inner regions 
such as azimuthal velocities close to that of light, large central redshift and ergoregions. The disk itself is 
likely to be stable against ring formation and the ratio of its binding energy to rest mass can be as large as 
1:10. The broad lines of how to construct all rotating disks arising from HKX-transforming the corresponding 
counter rotating model into a fully self-consistent model with proper rotation and partial pressure support 
has been sketched. It should be a simple matter to implement this method with the additional requirement 
that the pressure remains isotropic with a sensible polytropic index. One could then analyse the fate of a 
sequence of gaseous disks of increasing compactness and relate it to the formation and evolution of quasars 
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at high redshift. Alternatively, the inversion method described in section 4.1 yields a consistent description 
of the detailed dynamics for all disks in terms of stellar dynamics. 
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A. Equation of state for a relativistic 2D adiabatic flow 

The equation of state of a planar adiabatic flow is constructed here by relating the perfect fluid stress 
energy tensor of the flow to the most probable distribution function which maximises the Boltzmann entropy. 
Synge [21] gives an extensive derivation of the corresponding equation of state for a 3D flow. These properties 
are local characteristics of the flow. It is therefore assumed that all tensorial quantities introduced in this 
section are expressed in the local Vierbein frame. 

Consider a infinitesimal volume element d$ defined in the neighbourhood of a given event, and assume 
that the particles entering this volume element are subject to molecular chaos. The distribution function F 
of particles is defined so that ^(R, P) ddd 2 P gives the number of particles in the volume d centered on R 
with 2-momentum pointing to P within d 2 P. 

The most probable distribution function F* for these particles is then given by that which maximises 
Boltzmann entropy S subject to the constraints imposed by the conservation of the total energy momentum 
and by the conservation of the number of particles within that volume. This entropy reads in terms of the 
distribution F 

S = -dd f F\ogFd 2 P. (A.l) 



The stress energy tensor corresponds to the instantaneous flux density of energy momentum through the 
elementary volume d$ (here a surface) 

where the 1/\S\ factor accounts for the integration over energy-momentum space to be restricted to the 
pseudo-sphere P a P a = m 2 . Indeed the detailed energy-momentum conservation requires the integration to 
be carried along the volume element 

J 6 (P a P a - m 2 ) d 3 P = J 6 (P° 2 - £ 2 ) d 2 PdP° = ^ , (A3) 

where £ = \/m 2 + P 2 . It is assumed here that momentum space is locally flat. Similarly the numerical flux 
vector reads 

d 2 P 

F*P a -j^ . (AA) 

The conservation of energy-momentum of the volume element dd implies that the flux of energy momentum 
across that volume should be conserved; this flux reads 

T a/3 np dd = dd f F*P a d 2 P = Const. , (A5) 



where np is the unit time axis vector. Keeping the population number constant provides the last constraint 
on the possible variations of S 

<f)«n a d$ = dd j F*d 2 P = Const. . (A6) 



Varying Eq. (A A) subject to the constraints (A.b) (A.6), and putting 8S to zero leads to 

(logF* + l)8F* = a8F* + \ a P a 8F* , (A.7) 

where a and X a are Lagrangian multipliers corresponding resp. to (A.b) (A.6). These multipliers are 
independent of P a but in general will be a function of position. The distribution which extremises S 
therefore reads 

F* = Cexp(X a P a ) . (A8) 
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The 4 constants C and A a are in principle fixed by Eqs. (A. 5) and (A. 6). The requirement for F* to be 
Lorentz invariant - that is independent of the choice of normal n a - and self consistent, is met instead by 
demanding that F* obeys Eq. (A.2) and (AA) , namely 

C J P«exp(A^)^ = ^, (A9) 

C J P a P {i exp (A„P") ^ = T a ? , (A10) 
where <f>" and T" 13 satisfy in turn the covariant constraints 



dx a 



and (A. 11a) 

0. (A116) 



Equations (A.9)-(AAl) provide a set of thirteen equations to constraint the thirteen functions C, \ a , <f>", and T" 13 . 
Equations (A.9) and (A. 10) may be rearranged as 

€ = C— . T'^C^. «.4.12, 

where the auxiliary function $ is defined as $ = J exp(A„_P a ) d 2 P/\£\; $ is best evaluated using pseudo polar 
coordinates corresponding to the symmetry imposed by the energy momentum conservation P a P a — m 2 = 0; 
these are 

P 1 = m sinh x cos <f> , 

P 2 = m sinh x sin <f> , (A. 13) 

P° = im cosh x , 

where the coordinate x is chosen so that the normal n a lies along x = 0. In terms of these variables, d 2 P/\£\ 
then reads 

d 2 P/\£\ = dP 1 dP 2 /m cosh x = m sinh. xdxd<f> • (^-14) 

Moving to a temporary frame in which A is the time axis ( A = i\ = i(— A^A") 1 / 2 , X k = k = 1, 2), $ 
becomes 

y 00 _ 27r 

$ = 27rm / exp (—mA cosh x) sinh x^X = ex P ( — ^) • (A. 15) 

Jo 

From Eq. (A9)-(A10) and (A. 15), it follows that 



7^ = ^exp(-A) 



{A 2 + 3A + 3}^ + {A+l}^ 
A A 



(A. 16) 



^ = ^exp(-A){A+l}A«. (A17) 
On the other hand, the stress energy tensor of a perfect fluid is 

T a P = (e +p )U a U f) + Po 6 a P , (A18) 
while the numerical surface density of particles measured in the rest frame of the fluid £ is related to <f>" via 

£ = -€U a ■ (A19) 
Eq. (A. 16) has clearly the form of Eq. (A. 18) when identifying 

U a = \ a /\, (A. 20a) 

eo+Po =^0- exp (-A) { A 2 + 3A + 3} , (A206) 

Po = ^ eX p(-A){A + l} . (A20c) 
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Eliminating X,C between Eqs. (A.l7)-(A.20c) yields the sought after equation of state 

s = 2p + 



(A21) 



where po = E/A. This is the familiar gas law given that 1/A is the absolute temperature. The equation of 
state Eq. (A. 21) corresponds by construction to an isentropic flow. Indeed, using the conservation equations 
Eq. (A. 11a) dotted with U a , and Eq. (A. 116) together with U a U a = —1 yields after some algebra to the 
identity 



ds 



(A 2 + 3A + 3) 
A(A + 1) 



dU c 
dx° 



1 dE 

E ds 



(A. 22) 



where d/ ds = (dx a / dx a )(d / dx a ) is the covariant derivative following the stream lines. Writing the r.h.s of 
Eq. (A. 22) as an exact covariant derivative leads to the first integral 



ds 

Given that po = E/A, it follows 



1 + A 



2 log (A) + log (1 + A)- log (E) 



0. 



1 , P° \ P° 
1+ EJ E 6XP 



TTW^) h = Const - 



(A. 23) 



(A24) 



which in the low temperature limit gives E 2 po = Const. This corresponds to the correct adiabatic index 
7 = (2 + 2)/2. 



B. Distribution functions for counter-rotating relativistic disks 

For counter rotating disks in the absence of gravomagnetic forces, the inversion for distribution functions 
given in section 4.1 can be carried in close analogy to the Newtonian method given by Pichon & Lynden-Bell 
(1995) [20], and is sketched here. The relativistic flow is characterised by its stress energy tensor, 
which, in Vierbein components, reads 



(5.1) 



From the geodesic equation it follows that 

W, = -l => pW = ^[ £ *-h* 7u / 7H -g u ] 1/2 , (5.2) 

VTtt 

where the line element in the disk is taken to be 

ds 2 = -j tt dt 2 + 7RR dR 2 + 7<t>4> d<f> 2 . (5.3) 

Differentiating Eq. (5.2) implies 

dP (R) de , , dh 

— 77T = 7T7 , while dfW = — , (5.4) 

where the differentiation is done at constant R and h , and R and e respectively. Combinations of Eqs. (5.1) 
lead to 

A = 1H 1 ' 2 (T(«) - rW) - T( RR A = [ [ L(c,h)dhde _ ^ (B ^ 

A = = I I ff*£> h ] dhd < . (5.56) 
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Note that Eq. (5.5a) corresponds to Tolman's formula [23] which yields be integration the total gravitating 
mass of the disk. Introducing 

TZ 2 = lH / ltt , e = l - (e 2 - 1) , * = | (1 - la) and /* (e, h) = , (5.6) 
gives, for Eqs. (5.5): 

f f f+ (e, h) dhde , f f h 2 f+(e,h) dhde . „ . 

A = / / i±^J and A = / / — . 5.7 

Eq. (5.7) is formally identical to Eq.(2.9) of Pichon & Lynden-Bell [20], replacing (A, A, H, e, <J) by 
(RTi, R 3 p^, f,R,e,ip). In fact (A, A, 7?., e, \?) tends to i? 3 p^, f,R,e,ip) in the classical regime. 

Introducing the intermediate functions F(h,lZ), which is chosen so that its moments satisfy Eqs. (5.5), and 
following the classical prescription of Pichon & Lynden-Bell [20] yields 

1 K r T (dP/dn) dn , °? (dP/dn) dn 

f(e,h) = - / / /ft / , V Jh =, (5.8) 

where lZ p (h,e), and 1Z a {h,e) are respectively the apogee and perigee of the star with invariants (h,e), and 
lZ e (h) is the inner radius of a star on a "parabolic" (zero e energy) orbit with momentum h. Eq. (5.8) 
provides direct and systematic means to construct families of distribution function for counter-rotating disks 
characterised by their stress-energy tensor in the plane. Note that in contrast the dragging of inertial frames 
requires to fix three moments of the velocity distribution to account for the given energy density, pressure 
law and rotation law. 



C. Distribution functions for rotating disks: Lorentzian in momentum 

A possible choice for F is a Lorenzian parametrised in width, p2, mean, pi and amplitude, po such that 
F (v r p,R)dv r f, = PoP2ir where p = anc | v$ = e v {x - u) / R c ■ (C.l) 

pi + ip-pi) y l ~ v l 

The explicit function of v$ is found using dp = dv$/(l — v"^) 3 / 2 . In this instance, Eq. (4.12) becomes 

T<«> e- = P -RP1 r d JL = po , {C .2a) 

K J -co lp 2 2 + (p - J 

T y 9 >e = / = G{Po,Pi,P2) , (C.2b) 

w yl + p 2 ip2 2 + (p-pi) J 

T ( H ) e 4v = ?°Pl r pI^P = Po(P2+P2 2 +P1 2 ) ^ 2 ^ 

* J-°° (1+p 2 ) 2 (p2 2 + (P ~ Pi) 2 ) (1+P2) 2 +Pl 2 

where G(po, Pi, P2) is given by Eq. (C.4) bellow. The inversion of Eqs. (C.2) for (po, pi, P2), together with 
Eqs. (C.l), (4.14) and Eq. (4.3c) yields the distribution function of a Lorenzian flow compatible with the 
energy distribution imposed by the metric (v, £, u>) and the cut z = f(R). This inversion is always achievable 
provided the energy conditions Eq. (3.12) are satisfied. 

More generally, families of distributions functions corresponding to azimuthal momentum distribution 
such as 2p 3 po it -1 j{p\ + (p — Pi) 2 ) 2 , 8p| po (37r) _1 /(p2 + (p — Pi) 2 ) 3 ■ ■ ■ are readily derived in a similar manner 
while differentiating the r.h.s. of Eq. (C.2) with respect to p\. These higher order distributions have all 
vanishing number of stars with azimuthal velocity close to the velocity of light. This parametrisation is 
illustrated for the Kerr disk on Fig. 7.5. 
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C.l. The quadrature for G(po, Pi, P2) 

Changing variables to u such that p = (u — l/u)/2 Eq. (C.2b) becomes 



G(po,pi,p2) = — — / o o , / o — : — o — ^ ( c - 3 ) 

t Jo 4p 2 m + (« - 1 - 2pi w) 



This integral is then carried and yields 

Po (P2 + ipi) log (-2ip2 + % Pi ~ 2 \/l + (-«>2 + Pi) 2 ] 
Q{po,Pi,P2) = , h 

2^1 + (-ip 2 +p 1 ) 2 7T 

Po {ip2 - Pi) log ^-2ip 2 + 2pi + 2 \/l + (-«>2 + Pi) 2 ^) Po(p2 - ipi) log ^2ip 2 + 2pi - 2 ^1 + («>2 + Pif^j 

2^/l + (-i P 2+Plf TT 2^/l + (ip 2 + Pl f 7T 

Po (ipi - P2) log (^2ip 2 + 2 Pi + 2 \/l + (1P2 + Pif^j 

2^/l + (ip 2 + Pl f 7T 

(C.4) 



D. The relativistic epicyclic frequency 

The equation for the radial oscillations follows from the Lagrangian 



C = -\Je 2v [i - c^J 2 - t'^-' 1 " (1 + f' 2 ) R 2 - e~ 2v <j> 2 R 2 , (DA) 

where ( ) stands here for derivation with respect to the proper time r for the star describing its orbit. In 
Eq. (DA), <f> and t are ignorable which leads to the invariants: 

^ = h 8 -± = -e. (D.2) 
The integral of the motion for the radial motion follows from Eqs. (D.2) and U^lJ ' ^ = — 1, namely: 

e 2 C -2, ^ + //2 ) R 2 + ^1 {h _ e(jj) 2 _ e 2 e -2 V = _ x 

This equation, together with its total derivative with respect to proper time provides the radial equation of 
motion, having solved for the angular momentum h and the energy e of circular orbits as a function of radius 
when equating both R and R to zero, the relativistic generalisation of the classical epicyclic frequency is 
defined here to be the frequency of the oscillator calling back linear departure from circular orbits. Hence 
the equation for radial perturbation reads 

6R + k 2 6R= 0, (D.3) 

which gives for k 2 : 

K 2 =e~ 2Z (Z" - 2Z' 2 ) - e 2 e~ 2L (L" - 2L' 2 ) - 

V 7 V 7 (D.4) 



e 2K [ (h - ew) 2 (K" + 2K' 2 ) - e (h - eu) (4A'V + u") + eV 2 ] 
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where ( )' stands in this Appendix for d/dR = d/dR + /' d/dz and L, K, andZ are given in terms of the 
potential £ and v by: 

L = C + lcWl + /' 2 , (D.5a) 
Z = ( + log + /' 2 - f , (I>.56) 
K = 2v - C, - log ^1 + /' 2 - log i?. (D.5c) 



Here e and h are known functions of R given by the roots of: 



R =0 => e-' 2L e 2 =(h-euj) 2 e 2N + e~ 2Z , (D.6a) 
R = =>• e" 2Z Z' = e-^L'e 2 ^- (D.66) 
(/ l -ew)e 2iv [(/i-ewjAT'-ew']. 



E. Correction to Lemos' Self-Similar Disk 

The result given in section 6.1 agrees with that of Lemos, having made the following agreed corrections to 
his solution: 

• Eq.(_L 3.14) should read: (k + n) cot ±(k + n)7r = —Airp^rr with a minus sign, because the integration of 
6 is carried downwards. This gives n + k > 1 for the disks with positive pressures. 

• Lemos gives the stress energy tensor integrated through the plane using coordinate increment rather 
than proper length. Our result gives the latter which has an extra factor exp(C — v) and gives the 
physical energy tensor per unit area. 

• The exponent of cos[(n + k)ir/2] is 2n(k — n)/(n + k) 2 rather than An(l — 2n)/(k + n) 2 . 



